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ABSTRACT 

We  analytically  study  the  influence  of  convection  caused  by  horizontal  heat  transfer 
through  the  sides  of  a  vertical  Bridgman  apparatus.  We  consider  the  case  when  the  heat 
transfer  across  the  side  walls  is  small  so  that  the  resulting  interfacial  deformation  and  fluid 
velocities  are  also  small.  This  allows  us  to  linearize  the  Navier-Stokes  equations  and  express 
the  interfacial  conditions  about  a  planar  interface  through  a  Taylor  expansion.  Using  a  no 
tangential  stress  conditions  on  the  side  walls,  asymptotic  expressions  for  both  the  interfacial 
slope  and  radial  segregation  at  the  crystal-melt  interface  are  obtained  in  closed  form  in 
the  limit  of  large  thermal  Rayleigh  number.  It  is  suggested  that  these  can  be  reduced  by 
appropriately  controlling  a  specific  heat  transfer  property  at  the  edge  of  the  insulation  zone 
in  the  solid  side. 
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1  Introduction 

In  a  vertical  Bridgman  apparatus,  a  cylindrical  container  (Fig.  1)  containing  a  melt  of  a 
binary  alloy  mixture  is  translated  downwards  from  a  hot  to  a  cold  zone  so  as  to  cause  so¬ 
lidification.  The  uniformity  of  composition  of  the  resulting  crystal  and  the  relative  absence 
of  crystal  defects  are  desirable  features  for  technological  applications.  Ideally,  these  can  be 
achieved  if  convection  is  eliminated  and  the  crystal-melt  interface  is  planar.  However,  in 
practice,  this  is  difficult  to  ensure.  Considerations  of  constitutional  supercooling  (morpho¬ 
logical  instability)  and  the  need  to  avoid  transient  effects  due  to  container  ends  require  a 
relatively  large  temperature  gradient.  On  the  other  hand,  a  completely  one  dimensional 
imposed  temperature  gradient  that  would  occur  (Tiller  et  aF)  when  the  cylinder  sides  are 
insulated  would  require  unrealistic  large  temperature  differences  between  the  two  cylinder 
ends  given  that  the  length  of  the  cylinder  has  to  be  large  enough  to  avoid  transient  effects. 
Thus,  one  is  forced  to  a  configuration  where  significant  heat  flux  occurs  through  the  cylinder 
sides.  However,  when  this  happens,  fluid  next  to  the  cylinder  ends  is  hotter  than  the  fluid 
at  the  center  leading  to  convection  for  any  value  of  the  solutal  or  thermal  Rayleigh  num¬ 
ber.  The  review  by  Brown^  discusses  the  Bridgman  problem  in  great  detail.  Other  review 
papers^"®  deal  with  various  aspects  of  directional  solidification  in  general. 

There  are  many  papers  in  the  literature  that  address  the  problem  of  onset  of  convection 
in  finite  and  infinite  geometries  (see  references  6-15  and  references  there  in).  In  these  cases, 
the  equations  allow  for  a  basic  quiescent  state  (no  fluid  flow)  that  is  stable  upto  a  certain 
critical  Rayleigh  number.  Such  convection  can  be  termed  natural  convection,  as  opposed 
to  the  induced  convection  caused  by  heat  transfer  through  the  side  walls  of  a  Bridgman 
apparatus.  Since  this  paper  concerns  situations  which  in  the  absence  of  horizontal  heat 
transfer  is  a  thermally  stable  configuration,  natural  convection  is  not  relevant  in  this  context 
unlike  that  of  an  otherwise  thermally  or  solutally  unstable  arrangement  where  possibility 
exists  of  a  "resonant”  response  leading  to  vigorous  convection. 

There  are  many  papers  dealing  with  the  Bridgman  problem  itself  (see  Brown,  1988  for 
references).  Nonetheless,  the  first  set  of  fully  consistent  calculations  that  account  for  con¬ 
vection  in  the  fluid,  and  its  coupling  with  thermo-solutal  field,  heat  diffusion  in  the  solid 
and  a  non-planar  interface  appear  to  be  due  to  Chang  &  Brown^®  and  Adornato  &  Brown'^. 
By  assuming  a  quasi-steady  growth  when  the  transient  scale  is  small,  they  calculate  the 
nonlinear  steady  states  numerically.  Later  on,  these  calculations  were  further  extended  by 
Kim  L  Brown^*  by  including  effects  of  heat  transfer  through  an  ampoule  that  surrounds  the 
cylinder  in  an  experiment.  Unste2wly  transient  effects  have  also  been  studied  by  Brown  & 
Kim’®  and  these  calculations  are  in  good  agreement  with  experiment.  These  calculations 
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suggest  among  other  things  that  the  classic  one-dimensional  modelling  through  the  Scheil 
equation^®  is  an  oversimplification  since  it  assumes  complete  and  through  mixing  all  over 
the  fluid.  It  is  known  for  instance  that  between  the  diflFusive  limit,  where  convection  can  be 
totally  neglected  and  the  thoroughly  mixed  limit  for  which  Scheil’s  equation  applies,  there  is 
an  imperfect  mixing  zone  at  moderately  large  Rayleigh  number.  In  this  range,  radial  segre¬ 
gation  and  the  interface  deformation  are  larger  than  at  very  large  Rayleigh  number.  Indeed, 
it  was  pointed  out  that  reduction  of  gravity  by  a  factor  of  10^  —  10^  compared  to  earth 
may  be  detrimental  to  growing  a  crystal  of  uniform  composition.  Despite  such  advances  in 
theoretical  understanding,  it  is  difficult  from  numerical  results  alone  to  get  a  global  under¬ 
standing  of  the  parameter  dependences  since  there  are  so  many  of  them  (See  Table  1).  The 
trends  in  a  certain  subset  of  parameter  space  need  not  reflect  the  trend  in  other  ranges  of  the 
parameter  space.  Thus,  there  is  need  for  analytical  results,  which  is  likely  to  be  of  limited 
validity;  nonetheless  it  can  be  complimentary  to  numerical  calculations. 

The  only  analytical  work  that  we  are  aware  of  that  is  relevant  to  convection  in  the 
Bridgman  apparatus  is  due  to  Brattkus  &  Davis^^ .  They  analyzed  a  two  dimensional  model 
where  the  vertical  dimension  is  far  larger  than  the  horizontal  dimension  and  the  heat  flux 
through  the  side  walls  is  assumed  small.  Brattkus  &  Davis^^  specifically  concluded  there 
was  no  necessary  relation  between  radial  segregation  and  interfacial  shape,  a  hypothesis  put 
forward  by  Coriell  &  Sererka^^  based  on  diffusion  alone.  The  numerical  results  of  Brown  and 
his  coworkers,  on  the  other  hand,  suggests  a  strong  correlation  between  the  two.  However, 
the  Brattkus-Davis  analysis  ignores  the  insulation  zone  shown  in  Figure  1.  This  insulation 
zone  length  is  known  to  be  an  important  parameter  from  prior  numerical  work  (see  Chang 
&  Brown^®,  for  instance). 

Here  we  consider  steady  state^  inside  a  vertical  cylinder  (Figure  1)  between  0  =  0  and 
z  =  02,  where  02  is  assumed  a  constant.  As  in  Brattkus  h  Davis*’,  the  horizontal  heat  transfer 
is  assumed  small  enough  so  that  both  the  fluid  velocities  and  the  interfacial  deformation  are 
small.  Prior  numerical  work  (Chang  &  Brown’®  for  instance)  show  that  the  latter  assumption 
is  valid  even  for  relatively  intense  convection.  The  fluid  velocity  on  the  other  hand  has  to 
be  small  so  that  nonlinearities  in  Navier-Stokes  equation  can  be  ignored.  This  assumption  is 
clearly  unrealistic  in  many  experimental  situations;  however,  as  discussed  in  section  7,  there 
are  reasons  to  believe  that  the  results  on  radial  segregation  and  interfacial  deformation  will 
hold  in  part  of  the  nonlinear  regime  as  well. 

The  above  assumptions  allow  us  to  linearize  about  a  basic  one  dimensional  state  and 
Taylor  expand  the  interfacial  boundary  conditions  about  the  original  planar  boundary.  By 

^In  reality,  this  is  a  quasi-steady  state  (see  Brown  (1988));  however  for  purposes  of  this  paper,  we  make 
no  such  distinction. 
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modifying  the  no-slip  boundary  condition  at  the  cylinder  side  walls,  we  find  a  modal  rep¬ 
resentation  for  each  of  the  stream  function,  radial  temperature  and  concentration  gradient 
in  the  melt  that  decouples  each  modes  and  reduces  the  problem  to  a  set  of  finite  ordinary 
differential  equations  for  each  mode.  By  making  use  of  general  expressions  for  the  perturbed 
temperature  and  concentration  fields  in  the  solid,  we  find  that  the  interfacial  conditions  can 
be  expressed  in  terms  of  effective  boundary  conditions  on  the  melt  variables  at  the  original 
planar  interface.  With  the  effective  boundary  conditions  on  a  planar  interface,  we  solve  the 
linearized  melt  equations  and  find  explicit  expressions  for  interfacial  deformation  and  radial 
segregation  in  the  asymptotic  limit  of  large  thermal  Rayleigh  number  Rt-  Among  other 
things,  we  find  that  there  is  a  boundary  layer  near  the  interface  that  scales  as  where 

each  of  the  solutal,  thermal  and  fluid  velocity  field  change  rapidly.  The  radial  derivative  of 
the  temperature  field  and  the  radial  segregation  in  the  crystal,  each  scale  as  provided 

m  «  The  coefficient  of  Rj^^^  in  each  scale  decreases  exponentially  with  the  dis¬ 

tance  of  the  interface  from  the  end  of  the  insulation  zone  in  the  solid  side.  In  a  certain  range 
of  parameters,  consistent  with  many  experimental  conditions,  our  results  suggest  that  radial 
dependence  of  the  interface  shape  and  concentration  will  be  roughly  given  by  the  Bessel  func¬ 
tion  Jo(Air),  where  Aj  is  the  first  positive  root  of  Jo(An)  =  0  (i-e.  zero  of  J\).  Further,  we 
find  that  in  this  case,  the  Coriell-Sererka  (Coriell  &  Sererka^^)  hypothesis  of  proportionality 
between  interface  slope  and  radial  segregation  is  approximately  valid,  though  the  constant 
of  proportionality  is  different  from  what  these  authors  find  with  diffusion  only.  In  this  case, 
we  also  point  out  a  specific  condition  on  heat  transfer  near  the  solid  end  of  the  insulation 
zone,  which  when  satisfied,  will  result  in  minimal  interfacial  slope  and  radial  segregation  for 
large  Rt. 

2  Mathematical  Model 

In  a  Bridgman  apparatus  (Fig.  1),  a  binary  melt  is  contained  in  a  cylindrical  container 
of  radius  o  that  is  translated  downwards  with  constant  velocity  —UaZ,  where  i  is  a 
unit  vector  along  the  axis  of  the  cylinder  that  opposes  gravity,  as  shown  in  Figure  1.  The 
concentration  of  solute  (one  of  the  two  components  of  the  binary  mixture)  at  the  top  of  the 
cylinder  z  =  Z2  is  c  =  cj.  The  top  of  the  cylinder  is  maintained  at  temperature  T  =  Tj, 
while  the  bottom  of  the  cylinder  is  at  temperature  T  =  Ti,  which  is  significantly  smaller 
than  the  melting  temperature  To  of  a  planar  interface.  The  density  of  the  melt  is  assumed 
to  be  p  =  po  at  temperature  Tq  and  concentration  C2.  This  will  be  the  reference  density.  We 
scale  all  lengths  by  o,  all  time  scales  by  a/Ua,  all  mass  scales  by  poa^  and  all  temperatures 
by  T2  —  To,  leaving  us  with  non-dimensional  quantities  of  Table  1.  Concentration  will  remain 
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unsealed,  as  it  is  already  nondimensional.  Given  the  fair  amount  of  algebraic  manipulation, 
we  prefer  to  use  this  set  of  parameters  as  it  keeps  the  equations  looking  simpler.  Off  course, 
in  the  common  Engineering  literature,  it  is  common  to  express  these  quantities  in  terms  of 
Peclet  number  Pe,  Reynolds  number  Re,  Prandtl  number  Pr,  etc.  For  the  benefit  of  the 
readers,  we  have  prepared  Table  2  that  expresses  our  non-standard  parameters  in  terms  of 
the  more  well-known  dimensionless  numbers.  In  section  7,  we  also  discuss  our  results  using 
the  more  standard  notation  to  benefit  the  reader  who  is  more  interested  in  the  concrete 
results  than  in  th?  analysis. 

Using  the  nondimensional  variables  in  Table  1,  the  top  end  of  the  cylinder,  corresponding 
to  z  =  22,  is  maintained  at  nondimensional  temperature  T  =  1  and  the  bottom  end  is 
maintained  at  temperature  Ti  with  Ti  <  0  .  The  temperature  gradient  is  assumed 
strong  enough  to  avoid  constitutional  supercooling.  The  heat  is  allowed  to  flow  through 
the  sides  of  a  cylinder  for  22  >  2  >  2//  in  the  melt  zone  and  0  <  2  <  2/  in  the 
crystallized  zone.  The  heat  transfer  across  the  side  walls  is  such  that  the  fluid  velocities  are 
slow  and  the  interface  deformation  from  a  planar  interface  is  small.  The  precise  limitations 
placed  on  the  size  of  the  heat  transfer  by  this  assumption  will  be  examined  later  in  section 
6.  For  the  present,  it  suffices  to  assume  that  heat  transfer  is  sufficiently  small.  We  introduce 
a  cylindrical  coordinate  system  (r,  2)  where  r  is  the  radial  direction.  No  azimuthal 
angular  variable  is  necessary  as  the  flow  variables  are  assumed  to  be  axi-symmetric.  The 


mathematical  equations  for  the  steady  solution  in  the  melt  are 

v-VT  =  K  V^T,  (1) 

u-Vc  =  D  c,  (2) 

u  •  V  u  =  —  Vp  -h  {aT  ^  {c  —  C2)}  g  z  1/  V*  u,  (3) 

V  •  v  =  0,  (4) 


where  T  denotes  the  nondimensional  temperature,  c  the  concentration  of  solute  (one  of 
the  components)  relative  to  the  total  (measured  in  molar  fraction)  and  v  is  the  relative 
melt  velocity  (see  Table  1)  and  v  denotes  the  nondimensional  kinematic  viscosity  (i.e 
inverse  Reynolds  number  based  on  Ua).  In  equations  (3)  and  (4),  we  have  invoked  the  usual 
Boussinesq  approximation,  where  the  density  variation  due  to  change  in  temperature  and 
concentration  from  reference  values  are  included  only  in  the  forcing  term  on  the  right  hand 
side  of  (3).  Here,  a  (  >  0)  and  /?  are  the  non-dimensional  coeflBcients  of  volumetric 
expansion  due  to  increased  T  and  c  respectively.  When  the  solute  density  is  smaller  than 
the  alloy,  >  0.  However,  there  is  no  restriction  on  the  sign  of  ^  in  our  current  analysis. 
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On  the  crystal  side,  since  there  is  no  fluid  motion  and  so 


V  = 


z. 


Thus  the  equations  for  temperature  and  concentration  fields  are  given  by 

dT, 


dz 

d  c, 
dz 


=  D.  v*  c,. 


(5) 


(6) 


(7) 


Now  comes  the  boundary  conditions.  Denoting  the  r  and  z  components  of  v  by  (u,  xv) , 
we  take  the  boundary  condition  on  fluid  velocity  components  at  the  sidewalls  r  =  1  as 


u  =  0, 
dw 


(8) 

(9) 


Eiquation  (9)  is  a  no-stress  condition.  A  no  slip  condition  to  =  —  1  would  be  more  realistic; 
however,  we  could  find  no  simple  basis  representation  of  the  solution  in  this  ciise.  We  suspect 
that  aside  from  the  changing  the  nature  of  a  boundary  layer  near  r  =  1,  the  change  of  this 
boundary  condition  will  be  have  no  global  effect  on  crystal  shape  and  radial  segregation,  at 
least  in  the  li'  it  of  large  Rayleigh  number. 

The  condition  of  no  mass  flux  through  the  side  walls  imply 

dc 


dr 


=  0. 


Also,  the  condition  of  heat  flux  imply 


^  =  0  for  /(I)  <  z  <  z//, 

€  qi  T  =  cq2  for  z//  <  z  <  Z2. 


On  the  solid  side,  at  r  =  1  , 


dT, 

dr 


=  0  for  /(r)  >  z  >  zi, 


while 


dr 

+  e  r.  =  <  qu  for  zi  >  z  >  0. 


(10) 

(lla) 

(116) 

(12) 

(13a) 

(136) 


Each  of  cqi  ,  €qi  ,  eg,,  and  will  be  taken  as  constants  that  characterize  the  heat 

transfer  through  the  side  walls.  In  the  common  engineering  literature,  egi  and  cq\,  are 
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Biot  numbers  related  to  heat  transfer  through  the  side  walls  in  the  melt  and  solid  region 
respectively.  For  instance,  in  equation  (16)  of  the  Adornato  &  Brown’^  paper,  the  notation 
Bi{z)  with  a  specific  choice  of  a  piecewise  constant  profile,  is  related  to  the  constants  f<7i, 
and  eqi,  while  Bi{z)6oo{2)  is  related  to  eq2,  and  eq2,,  provided  we  ignore  the  ampoule  in 
their  model.  The  quantity  f  will  be  a.ssumed  small,  while  each  of  qi  and  92  assumed  0(1). 
e  will  be  our  perturbation  parameter.  The  precise  multiplicative  decomposition  of  the  Biot 
number  egi  into  q\  and  e  is  unimportant  since  in  the  final  results,  only  the  product  eq^ 
appears.  The  same  is  true  for  eq2  ,  eq^^  and  eq2,- 

At  the  cylinder  top,  i.e.  z  =  22  >  we  have  the  nondimensional  temperature 


T  =  1,  (14) 

C  =  C2,  (15) 

u  =  0,  (16) 

w  =  -1,  (17) 

with  C2  being  assumed  a  constant.  Similarly  at  the  bottom,  z  =  0  ,  we  get 

T,  =  r,,  (18) 

c,  =  C2,  (19) 


where  <  0  and  assumed  smaller  than  the  nondimensional  melting  temperature  at  the 
interface.  Note  that  one  can  use  more  general  boundary  conditions  than  (19),  however  this 
appears  to  be  most  relevant  since  in  the  limit  of  an  infinitely  long  cylinder,  conservation  of 
mass  dictates  that  the  concentration  should  be  the  same  at  z  =  ±00  . 

Now  comes  the  boundary  conditions  on  the  solid  melt  interface  z  =  f{r)  .  We  will 
assume  local  thermodynamic  equilibrium  since  the  relaxation  time  of  the  departure  from 
equilibrium  can  safely  be  assumed  to  be  much  smaller  than  the  typical  time  scales  in  this 
problem.  Thus,  at  the  interface,  we  satisfy  the  melting  condition  given  by  the  solidus- 
liquidus  line  with  the  incorporation  of  the  Gibbs-Thompson  effect  (i.e.  lowering  of  the 
melting  temperature  in  the  presence  of  curvature  effects) 

{/"  1  f  ] 

- nn  + - TTH  f  •  (20) 

(1  +  r  (1  + 

Further,  from  the  solidus  liquidus  line,  we  must  have 

c,  =  k  c.  (21) 
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(22) 


The  continuity  of  heat  flow  implies  that 

.  dT 
I  V  •  n  —  K  — -  = 
on 


■Tg  z  n  —  Kg 


dn 


+  St  V  ■  n, 


where  n  denotes  a  unit  vector  normal  to  the  melt-solid  interface  pointing  towards  the  liquid, 
and  ^  denotes  component  of  gradient  in  that  direction,  St  is  the  nondimensional  latent  heat 
(Stefan  number).  From  the  conservation  of  solute  concentration  across  the  interface,  we  must 
have 


c  V  -  n  -  D  —  -Cg  z 
an 


n 


Dg 


dcg 

dn 


(23) 


Further,  from  integrating  (4)  over  a  small  control  volume  centered  about  a  point  on  the  in¬ 
terface,  it  follows  that  normal  velocity  is  continuous  when  the  difference  of  densities  between 
the  melt  and  the  crystal  is  neglected.  Thus, 

V  •  h  =  —  z  ■  h.  (24) 


Difference  of  densities  between  solid  and  liquid  side  of  the  interface  can  be  accommodated 
by  including  an  additional  —  •  n  term  on  the  right  hand  side  of  (24);  however,  here 

we  will  be  primarily  interested  in  situations  where  this  is  not  important.  Further,  at  the 
interface,  if  we  assume  a  no-slip  condition  for  the  fluid  flow  relative  to  the  solid,  then 

u  -f  (tn  -h  1)/'  =  0.  (25) 

While  the  no-slip  condition  (25)  at  the  melt-solid  interface  seems  to  be  widely  used  in 
the  literature,  we  are  unaware  of  any  convincing  physical  argument  this  is  superior  to  a 
no  tangential  stress  boundary  condition.  It  may  be  pointed  out  that  unlike  the  case  of 
a  fluid  next  to  a  solid  where  good  experimental  evidence  exists  for  the  no-slip  boundary 
condition,  the  liquid  molecules  in  this  case  do  not  preserve  their  identity  as  it  goes  through  a 
phase  transition.  Nonetheless,  despite  this  uncertainty,  we  find  that  the  eisymptotic  scalings 
and  parameter  dependence  presented  later  in  this  paper  remain  the  same  with  a  no-stress 
boundary  condition,  only  the  scaling  constants  differ. 

Note  that  the  mathematical  model  with  equations  (l)-(4),  (6),  (7)  and  boundary  and 
interfacial  conditions  (8)-(25)  contain  twenty  irreducible  nondimensional  parameters  Z),  Dj, 
«,  Kg,  ag,  pg,  v,  Z2,  zj,  zjj,  tq\,  tqi,,  tq2,  tq2,,  C2,  T\,  do,  m,  k  and  St-  The  relation  of  the 
less  standard  parameters  in  this  list  with  the  conventional  Peclet  number,  Prandtl  number, 
Reynolds  number,  etc.  is  given  in  Table  2. 


3  Steady  state  for  zero  epsilon 

Despite  the  complexity  of  the  general  equations,  there  exists  well  known  simple  solutions  to 
the  above  equations,  as  determined  originally  by  Tiller  et  aP.  If  e  =  0  ,  a  quiescent  state 
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is  a  simple  solution  to  the  above  equations  in  which  the  melt  velocity  v  =  —  z  and  the 
melt-solid  interface  is  planar,  i.e.  /  =  zq  ,  &  constant.  The  temperature  and  concentration 
fields  in  this  case  are  denoted  by  a  superscript  0  as  they  are  the  leading  term  of  an  expansion 
for  small  e  .  They  are  given  by 


On  the  solid  side,  zq  >  z  >  0  , 


-  l) 

1  -  Ti 

(26) 

V 

_  i’ 

C2  -  Ci 

(27) 

1 

7 

1 

1 

_ 

i) 

^  -  l’ 

(28) 

C2  -  k  a 

f  -  r 

(29) 

where  c,-  is  the  concentration  value  on  the  melt  side  of  the  interface  and  Ti  is  the  interfacial 
temperature.  To  determine  the  three  constants  Ti  ,  c,  and  zq  ,  we  use  the  boundary 
conditions  (20),  (22)  and  (23)  which  in  this  case  simplifies  to 


Ci  - 


Ti  = 

-  m  c<, 

(30) 

1  -  Ti 

- 

(31) 

_  1 

+  Of, 

-  1 

C2  -  Ci 

s=  k  Ci 

C2  —  k  Ci 

(32) 

_  1 

-  1 

In  this  case,  the  pressure  p  =  p^(z)  is  hydrostatic  and  satisfies 

^  =  ojr”  +  ^(,(c»-c,). 


(33) 


Note  that  all  other  equations  and  boundary  conditions  are  trivially  satisfied  by  the  solutions 
(26)-(29)  and  (33)  provided  they  satisfy  (30)-(32). 

An  important  limiting  case  (valid  in  many  experiments)  is  that /c/22  >>  li«s/2;2  >>  1- 
This  simplifies  expression  for  T®  and  7^  as 


T°  ~  -me,  -f  — — ^(1  +  me,), 

Zj  -  Zo 

~  -me,  -  — — —{Ti  +  mci), 
^0 


(34) 

(35) 
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each  of  which  is  a  linear  function  of  z.  In  addition,  when  (zj  —  ^o)  >>  D,  St  <<  n/z^, 
St  <<  /c,/22  and  Da  <<  20)  (conditions  that  appear  to  be  valid  in  some  experiments), 
the  matching  conditions  imply  that  the  interfacial  concentration 


c.  ~ 

(3G) 

.  f, ,  «(i  +  T‘;2)  r' 

'1  ■ 

(37) 

Notice  that  the  quantity  is  the  ratio  of  the  temperature  difference  between  the 

top  and  the  interface  and  the  temperature  difference  between  the  interface  and  the  bottom. 
Clearly,  by  controlling  the  location  of  the  insulation  zone,  zq  can  be  made  to  lie  between 
z  =  zj  and  z  =  2//,  as  will  be  assumed  here. 

Until  we  get  to  a  discussion  of  the  concrete  formulae  (212)  and  (213)  for  interfacial  shape 
and  radial  segregation  in  sections  6  and  7,  the  only  simplification  that  will  be  used  is  that 
each  of  and  is  a  linear  linear  functions  of  2  as  in  (34)  and  (35),  since  /c/22  and  Ka/z2 
are  assumed  large.  Otherwise,  the  analysis  will  proceed  with  the  assumption  that  each  of 
c“,  T®,  c“  and  Tf  are  known  from  (26)-(32). 


4  Perturbed  Steady  state  for  nonzero  e 

Now  consider  a  small  nonzero  c  .  In  this  case,  the  dependent  variables  can  no  longer  just 
depend  on  2  .  As  is  well  known,  the  presence  of  a  radial  thermal  gradient  means  that  the 
quiescent  state  is  no  longer  a  steady  state  solution  to  the  problem.  We  express  solutions  as 
a  perturbation  expansion  in  powers  of  e  : 


r  =  r®  +  c  r’  +  .., 

(38) 

c  =  c®  4-  C  C*  +  .., 

(39) 

«  =  c 

(40) 

w  =  — 1  4-  c  «7*  4- 

(41) 

p  =  p®  4-  ep'  4- 

(42) 

Ta  =  +  .., 

(43) 

Cs  =  c°  +  ec]  +  .., 

(44) 

/  =  zo  4-  e  4-  —  . 

(45) 
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Substituting  these  into  (1-4),  (6)  and  (7),  and  equating  the  t  terms  in  the  resulting  equations, 


-  ^  -1-  «)*  r®'  =  «  T*, 

dz 

(46) 

-^  +  „'c»'  =  DV«c', 
dz 

(47) 

— 

du^  dp^  ^21  i 

dz  dr  r^ 

(48) 

dw^ 

dz 

-  ^  -1-  (  O!  +  )9  c*)  g  -h  V  w^, 

dz 

(49) 

;  |(-’)  +  I-’  =  »■ 

(50) 

On  the  solid  side,  we  have 

dT}  „2 

az  = 

(51) 

=  D.  c). 

dz 

(52) 

We  note  that  in  this  coordinate  system,  the  operator  is  given  by 

2  _  ^  1  ^  ^ 

^  r  dr  dz^ 

The  boundary  conditions  on  r  =  1  for  zq  <  ^  <  ^2  are 

u*  =  0,  (53) 

^  =  0.  (54) 

%  =  0,  (55) 

dT' 

=  0  for  zo  <  z  <  zji,  (56a) 

or 

dT^ 

T®  =  02  for  zji  <  z  <  zj.  (566) 

dr 

On  r  =  1  for  0  <  z  <  Zo  ,  from  (12)-(13),  we  have  the  boundary  conditions 

dc^ 

ir  = 

dT^ 

=0  for  Zo  >  z  >  ^/,  (58a) 

dr 

dT^ 

^  -  91.  +  92.  for  zi  >  z  >  0.  (585) 
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The  boundary  conditions  at  the  top  2  =  23  are 


II 

0 

(59) 

c’  =  0, 

(60) 

u*  =  0, 

(61) 

=  0, 

(62) 

and  at  the  bottom,  2  =  0, 

r;  =  0, 

(63) 

c]  =  0. 

(64) 

On  the  original  interface  2  =  20  ,  (20)  impli^ 

r;  +  fT^’  =  T'  +'  /*  T®'  =  -mc°'  f  -  mc^  +  do  + 

(65) 

From  (21)-(25),  we  get 

c°'  =  k  c'  +  k  c°', 

(66) 

(67) 

-D  c”"  /■  -  0  ^  +  (t  - 1)  [c‘  +  c""  /'  ]  =  -D,  cf  /'  - 

(68) 

II 

0 

(69) 

=  0, 

(70) 

respectively.  We  can  eliminate  pressure  from  (48)  and  (49).  Also,  from  (50),  it  follows  that 
=  (u*,  ti;^)  =  curl  A  for  some  vector  field  A  that  can  be  chosen  to  be  divergence  free.  For 
the  axisymmetric  flow  under  consideration  A  can  be  reduced  to  a  scalar  ‘stream  function’ 
xj)  so  that  the  velocity  components  of  can  be  written  as 


It  =  — 


w'  =  - 


1  d{rxp) 
dr 


Then,  from  (47)  and  (48),  we  get 


=  v(c  + 


a* 

dz^ 


\  *  ar*  dc^ 

)  ^  +  + 


dr 


(71) 

(72) 


(73) 
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where  the  differential  operator  L  is  defined  such  that 


£  u  = 


fl  Q{ru)  1 

\  r  dr  ) 


On  taking  the  partial  derivative  of  each  of  equations  (46)  and  (47)  with  respect  to  r  ,  we 
obtain 

On  taking  partial  derivative  of  (51)  and  (52)  with  respect  to  r,  we  get 


-LEI  =  ^  /'/•  +  ^ 'I  ^ 

dz  dr  *  \  dz'^  )  dr  ' 

dz  dr  ~  ‘  \  ^  dz^  )  dr' 


The  boundary  conditions  (53)  and  (54)  imply  that  without  any  loss  of  generality, 

V»  (1,^)  =  0, 

£  rlf{l,z)  =  0, 

while  the  boundary  conditions  (61), (62),  (69)  and  (70)  imply  that 


zj)  =  0, 

(81) 

^  = »' 

(82) 

V'  (rjZo)  =  0, 

(83) 

II 

p 

(84) 

Further  on  taking  the  tangential  derivative  (i.e.  derivative  with  respect  to 
equations  (65)-(68),  we  obtain 

r  )  of  each  of  the 

,  rV  rpO>  1  rl'  rpO> 

ar  ’ 

(85) 

^  +  /-T”'  =  -me”'/"  -m^+doCf<', 

(86) 

(87) 
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(88) 


-« r”"  /■'  -  » I-  — - ''' 


dz  dr 


=  -K,  r 


d  dT] 


dz  dr  ’ 


— +  c®' 
dr  ^ 


=  -B.  cr  /■'  -  (89) 


In  addition  to  the  above  boundary  conditions,  elementary  considerations  of  smoothness  of 
each  of  the  variables  T*  ,  ,  c*  ,  c]  and  v*  in  the  neighborhood  of  r  =  0  together 

with  consideration  of  axisymmetrical  flow,  leads  us  to  conclude  that  as  r  — »  0  , 


v»  ~  <1  (z)  r. 

(90) 

dT^ 

dr  ~ 

(91) 

dT^ 

^  ~  H.)  r. 

(92) 

d(^ 

*  ~ 

(93) 

^  ~  ‘sW’-. 

(94) 

o,.  constant  r 

(95) 

for  some  functions  U  through  U  >  their  precise  form  being  unimportant, 

5  Series  representation  of  solution 

The  problem  (equations  (73),  (75)-(78))  in  conjunction  with  (55)-(60),  (63),  (64),  (79)-(89) 
and  assumptions  (90)-(95)  completely  determine  each  of  the  unknown  functions  , 

,  -^c]  and  V’  ^  functions  of  r  and  z  and  p'  as  a  function  of  r. 

We  expand  each  of  the  variables 


dT^ 

dr 

oo  2 

”  S  JiM 

(-2^) 

r), 

(96) 

dc' 

dr 

a  ~  2 

■  0  k  MK) 

bn  {z)  Jl(A, 

^  r), 

(97) 

ij)  = 

90^  2 
k  JiM 

Cn  (z)  «/l(An 

r), 

(98) 

dTl 

dr 

OO  2 

“  S  JUK) 

On  (2)  Jl{K 

r), 

(99) 

dr 

a  ^  2 

bn  {z)  Ji(A, 

»  r), 

(100) 
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oo  o 

/’’  =  E  -pm 

n=l 

where  A„  is  the  positive  zero  of  the  Bessel  function  Ji  .  Notice  that  this  is  a  suitable 
representation  since 

C  MXn  r)  =  -A^  J,  (A„  r) 

and  the  boundary  conditions  (55),  (57),  (79)  and  (80)  are  automatically  satisfied  by  each 
mode.  The  representation  of  ^  and  ^  in  (96)  and  (99)  may  appear  to  contradict  (56b) 
and  (58b).  This  is  not  the  case  since  the  convergence  of  the  series  in  (96)  and  (99)  for 
z  >  zii  and  z  <  z/  is  only  in  the  mean  and  not  pointwise  at  r  =  1.  In  this  range  of  z, 
there  is  a  slow  ^  decay  of  and  as  r»  — »  oo,  which  implies  that  one  cannot  calculate 
the  derivative  of  the  series  (96)  and  (99)  term  by  term.  Thus,  it  is  not  possible  to  plug  the 
series  expressions  directly  into  the  differential  equations  and  obtain  the  correct  equations 
for  o„  and  am-  Instead,  we  multiply  each  of  (73),  (75)-(78),  (85)-(89)  by  r  (A„  r)  and 
integrate  with  respect  to  r  from  0  to  1.  On  integrating  by  parts,  and  using  (55)-(58),  (79), 
(80),  (90)-(95)  and  the  following  relations: 

“”W  =  /o' 

'"W  =  77rji(K)l  ^rMKr)Hr,z), 

=  Jik)  /«’  *  ■■  ^ 

/o’  *■ "  ^ 

we  obtain  three  ordinary  differential  equations  for  each  n  on  the  melt  side  and  two  ordinary 


differential  equations  on  the  solid  side.  These  equations  can  be  written  as 

£2  K  =  Rc  h{z)  Cn  (oT  Zo  <  z  <  Z2  ,  (102) 

C3  Cn  =  -On  -  bn  foT  Zq  <  Z  <  Z2  ,  (103) 

On  =  -A^  Rt  Cn  +  A„  [  T®  +  92  ]  for  zn  <  z  <  Z2,  (104a) 

=  -\l  Rt  Cn  ioT  Zo  <  z  <  zii,  (1046) 
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where 


and  the  differential  operators 


*2,  =  0,  for  0  <  z  <  zo  , 

(105) 

=  [  -91.  +  92,  ]  for  0  <  z  <  zi 

(106a) 

jCi,  =  0  for  zi  <  z  <  zq 

(1066) 

1 

T 

II 

'TT 

(107) 

-  KU  ' 

(108) 

gfi  /  c,  -  C2  \ 

^  -  Dv  U  - 

(109) 

i  —  —  A*  C.  =  —  4.  —  —  - 

K  dz  “  ’  dz^  Kt  dz 

(110) 

(111) 

_  A'  +  l(£.  _  aA  a 

[dz^  V  ^  j/  [dz^  "j  dz’ 

(112) 

^  ^'dz2  dz  ” 

(113) 

Note  that  since  kJz^  »  1,  (34)  holds  and  therefore  is  a  constant.  So  the  thermal 
Rayleigh  number  Rr  defined  in  (108)  is  also  a  constant  (Note  the  definition  in  terms  of 
dimensional  variables  as  well  in  Table  2). 

From  (59)-(60),  (81)  and  (82),  we  find  that  at  the  cylinder  top, 


aniz2)  =  0, 

K  (^2)  =  0, 

Cn  (^2)  =  0, 

<  (22)  =  0. 

The  boundary  conditions  (63), (64)  at  the  cylinder  bottrnn  imply 

a«i(0)  =  0, 

(0)  =  0. 

At  z  =  zo  >  we  find  from  (85)-(r^  tha^ 

®»n(^^o)  +  i^o)  =  <ln(2o)  +  T^{zo), 


(114) 

(115) 

(116) 

(117) 

(118) 

(119) 

(120) 
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an(^o)  +  dn  T^izo)  =  -m  C^\zo)  da  -  m  ^  bn{zo)  -  do  Xl  da,  (121) 

b<m{zo)  +  da  -  c°\zo)  =  k  ba{zo)  +  k  da  -  c“'(«o),  (122) 

a  a 

-K  a!^  (zo)  -  K.  da  T^\zo)  =  -«»  a!^  (-Zo)  -  «.  da  Tf\zo),  (123) 

-Db'a{Zo)  -  Dda-C^\zo)-\-{^-\)[ba{Zo)+C^\zo)-da]  =  -D,b'^  {zo)- D,da-cf  {zo).  (124) 
a  a  a 

FVom  (83)  and  (84),  it  follows  that 

Ca(zo)  =  0,  (125) 

<  (zo)  =  0.  (126) 


Before  proceeding  further,  for  purposes  of  reducing  algebra,  we  simplify  the  equations 
further  by  assuming  that 

(a)  «,  k/z^  ,  and  >>  1 

(b)  do  =  0 

With  the  assumption  (a)  as  above,  we  replace  equations  (104a,b)  and  (106a,b)  by  simpler 
equations 

a„  =  -A^  Rt  Ca  +  \a[-qi1^  +  gz]  for  zjj  <  z  <  Z2  ,  (127a) 


£l  Oa  =  -A^  Rt  Ca  ioT  Zo  <  z  <  Zjl  ,  (1276) 

£i  =  A„  [  -9i,  +  92.  ]  for  0  <  z  <  zi  ,  (128a) 

C\  =  0  for  2/  <  z  <  ^o  »  (1286) 

where 

A  =  53  -  a;.  (129) 

The  interfacial  matching  conditions  (120)-(124)  couple  the  concentration  and  temperature 
fields  in  the  solid  side  to  the  melt  variables.  Using  (105),  (118),  (119)  and  (128),  expression 
for  Can  and  b^a  can  be  obtained  in  terms  of  two  arbitrary  constants.  By  eliminating  these 
arbitrary  constants  between  a,„(2;o),  a^n(^o)  and  bsn{zo),  6^(^o)j  the  matching  conditions 
(120)-(124)  can  be  written  as  two  effective  boundary  conditions  a.t  z  =  zq  on  the  melt 
variables.  As  shown  in  Appendix  A,  these  boundary  conditions  are 

-^<(^o)  +  Aan(^o)  +  ^zbnizo)  =  (130) 

b'^{zo)  +  miba{zo)  +  m2aa{zo)  =  0,  (131) 

where  the  effective  interfacial  parameters  ^2i  ”^2  and  /?3  are  defined  as: 


«,(r°'  +  mc»'(zo))(l  +  e-2^"^) 

«(mc“'(zo)  +  7’°')(1  -  e-2^"^)  ’ 


(132) 
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From  (102),  (103)  and  (127),  we  get  an  eight  order  linear  differential  equation  for  c,,  : 

[£i£2-C3  +  XlRcCih{z)  -  XlRTC2]cn  =  £2  [-A„[-9iT®  +  92]]  for  zn  <  z  <  Z2, 

(143a) 

[£i£2£3  +  X\RcC\h{z)  —  c„  =  0  for  zo  <  z  <  zn.  (1436) 

Alternatively,  instead  of  one  equation  for  c„,  we  can  write  the  following  equations  for  and 
On 

[£3£2  -  DX^Rt  +  XlRch{z)]  Cn  +  oj,  =  -£)An[-9ir®  +  92]  for  zn  <  z  <  22,  (144a) 

[£3£2  -  DXIRt  +  A^/2c^(-2)]  Cn  +  aJj  =  0  for2o<2<  zjj,  (1446) 

£i«n  =  -XIRtc^  +  An[-9ir®  +  92]  for  zn  <  z  <  22,  (145a) 
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(1456) 


C\an  =  —XlRrCn  for  zq  <  z  <  zji  . 

The  boundary  conditions  are  (125),  (126),  (130)  and  (131)  at  z  =  zq  and  (114)-(117)  at 
z  =  Z2.  In  addition,  we  require  that  a„  and  its  first  derivative,  6„  and  its  first  derivative 
and  c,i  and  all  its  first  three  derivatives  are  all  continuous  at  2  =  Z]j.  From  (103),  this  is 
equivalent  to  requiring  that  a,i  and  its  first  derivative,  c„  and  its  first  five  derivatives  are 
all  continuous  at  z  =  zjj.  These  are  enforced  on  a  suitable  representation  of  the  general 
solution,  as  will  be  discussed  shortly. 

Since  (143a)  is  an  eighth  order  non  homogeneous  linear  equation,  we  can  express  the 
general  solution  for  zjj  <  z  <  2j  in  the  form 

Cn  =  c^(^)  +  't  + i;  (146) 

j=l  j=3 

where  Cnp{z)  is  a  particular  solution  to  (143a),  Fj  are  eight  independent  solutions  to  the 
homogeneous  equation  (143b)  and 

7=  (147) 

Factors  containing  powers  of  7  in  (146)  can  be  absorbed  as  part  of  arbitrary  constant  Cj 
by  suitable  redefinition.  We  choose  not  to  do  so  for  we  want  to  explicitly  show  the  Rj  scaling 
of  the  eventual  answer.  With  our  choice  of  it  will  turn  out  Cj  =  0{\)  as  Rr  — >  00. 
Questions  of  determination  of  functions  F/  and  Onp  will  be  set  aside  for  the  moment.  Once 
Cn  is  obtained,  expression  for  a„  and  6n  be  found,  at  least  in  principle  from  (144),  (145)  and 
(103).  We  write  these  symbolically  in  the  form 

an  =  ar^{z)  +  53  l~^CjF^{z)  ,  (148) 

3  A*  /L  * 

*«  =  ^,  +  E  7-'c,pf(z) + -^  E  .  (1«) 

>=1  '  J=3 

where  F“,  F*  are  related  to  F/;  while  particular  solutions  a„p  and  6„p  are  determined  in 
terms  of 

For  Zo  <  z  <  zji,  i.e.  in  the  insulated  zone,  we  write 


3  8 


'Er%Fi +'£r^CiF; 

i=i  j=3 

(150) 

8 

j=i 

(151) 
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(152) 


=  1^7 

i=i 


-1 


XI R, 


j=3 


In  appendix  B,  we  derive  asymptotic  representation  for  each  of  FJ,  F“,  Fj,  Cn^,  a„^ 
and  bnp  for  large  Rt  that  is  uniformly  valid  for  all  An,  provided  |/2c|  <<  Rt-  We  also 
derive  simpler  expressions  for  the  special  case  Rf  »  1  with  7®  >>  A*/lc  and  this 
latter  expressions  will  be  quoted  here  as  these  are  the  only  ones  needed  in  obtaining  concrete 
results  for  l^cl  <<  Fj*- 

In  the  special  case,  A^|Fc|  «  7®,  the  leading  order  asymptotic  behavior  is  given  by: 


~  (ajat  +  as)  4a;at  +  a‘P’ 

pc  ^  £1  pPi(*-*o)  pc  ^  ?1  pPa(*-*2) 

F;  ~  for  3,4,5  , 

F/  ~  for  j  =  6,7,8  , 

where  pi  and  pa  are  given  by 

~W  ^ 

Also,  7j  for  j  =3  to  8  are  the  six  independent  roots  of 

-  aJAt  =  0, 


(153) 

(154) 

(155) 

(156) 

(157) 

(158) 


where 

i.b)  =  »“  -  Ai, 

t3(sl  =  (»"  -  Ail)“  +  -  A^). 

The  roots  of  (158)  are  labelled  such  that  for  Rt  »  AJ[, 

7j  ~  7‘*^i  (159) 


where 


U3  =  =  -078 

U4  =  — 1  =  — Wy 
UJ,  =  -e-'*'’  =  -U)6 
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(160). 

(161) 

(162) 


o„p  is  general  is  a  complicated  expression,  given  in  the  appendix.  We  will  not  need  the 
general  expression  in  what  follows  except  to  not  that  in  the  special  case  when  <<  Rj, 
the  general  expression  simplifies  such  that 

~  XlRr  +  AS  ^[XIRt  +  Xt¥  LA-D-'] 

where 

L2[y\  =  Diy^  -  \l)  +  y. 

For  7®  >>  A^|/Zc|  (A„  unrestricted  otherwise), 


(164) 

- J=  3,4,5  , 

(165) 

i  =  6,7,8  . 

(166) 

The  asymptotic  form  of  6„p  is  simply  obtained  by  using 

^np  —  ®np  ,^3[^p] 

(167) 

For  7®  » 

A^|/2c|,  The  independent  homogeneous  solutions  for  6„  have  the  asymptotic 

behavior 

pb  ^  _  gPl(z-Zo)  ^  pb - gP3(*-*2)  ^ 

(168) 

F?  ~  -h  ..  for  i  =  3,4,5  , 

(169) 

^  +  „  for  >  =  6,7,8  , 

(170) 

where 

1 

7i=7i-;p 

These  asymptotic  results  for  a  special  case  will  be  used  later  in  obtaining  concrete  ex¬ 
pressions  for  quantities  of  interest.  It  is  to  be  noted  that  the  asymptotic  expressions  (155), 
(156),  (165),  (166),  (169),  (170)  become  invalid  when  |A*/?c|  is  the  same  order  as  7®  or  larger. 
This  can  occur  even  for  Rc  «  provided  A„  is  large  enough.  In  Appendix  B,  section 
2,  we  derive  more  general  expressions  for  the  solu'^’on  of  the  homogeneous  equation  will  be 
uniformly  valid  for  all  A„  provided  Rt  «  Rc.  The  expressions  (155),  (156),  (165),  (166), 
(169)  and  (170)  will  have  to  be  replaced  by  these  expressions.  The  remaining  asymptotic 
relations  (153),  (154),  (163),  (164),  (167)  and  (168),  however,  remain  valid  in  this  range  of 


np 


L2[-D-^]\ 


(163) 


parameters. 
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Despite  differing  asymptotics  in  different  regimes,  at  this  stage  we  prefer  to  think  of 
relations  (146)-(152)  as  exact.  The  only  eispects  of  the  asymptotic  results  of  Appendix  B  in 
differing  regimes  that  will  be  used  in  this  part  of  our  analysis  is  that  in  all  cases, 

(i)  Each  of  and  have  a  common  expor«ntial  part  e**^^*^  multiplied  by  an  algebraic 

dependent  term  that  depend  on  superscript  a,  b  or  c.  For  instance,  when  <<  7^, 


=Pi{z-zq),  W2  =  P2{z-Z2),  Wj  =  7j(2-zo)  for  j  =  3,4,5  Wj  =  7j(z-22)  for  j  =  6,7,8 

(171) 

(ii)  The  variation  in  z  of  for  j  =  3, 4. .8  is  much  larger  than  the  variation  of  the  algebraic 
prefactors  when  Rj  »  1,  while  for  A„  >>  1,  the  variation  of  the  exponential  part 

for  j  =  1,2  dominate  the  variation  of  the  prefactors. 

(iii)  For  j  =  1 , 3, 4, 5,  Re  Wj  >  0,  while  for  j  =  2, 6, 7, 8,  Re  Wj  <  0. 

(iv)  The  choice  of  the  constants  Cj  and  Cj  and  the  prefactors  is  such  that  Wj{zQ)  —  0  for 
j  =  1,3,4, 5,  while  Wj{z2)  =  0  for  j  =  2,6, 7,8. 

As  a  consequence  of  the  above  properties,  it  is  clear  that  for  large  /?r,  we  have  each  of 
Fj{zQ),  Ff{zQ)  and  Fj(zo)  will  be  exponentially* small  in  Rt  for  j  =  6,7,8,  while  Fj{z2), 
F^{z2)  and  (22)  are  exponentially  small  in  Rt  for  j  =  3,4, 5.  Again  for  An  >>  1,  each  of 
F2“(2o),  ^2(20)  ,  ^2^(20),  F“{z2),  F*(22)  and  F^{z2)  is  exponentially  small  in  An.  These  facts 
will  be  needed  later  in  simplifying  matrix  equations. 

The  continuity  of  Cn  and  its  first  derivative,  Cn  and  its  first  five  derivatives  can  be  written 
as  the  following  matrix  relation 

GX  =  R,  (172) 


where  G  is  a  8  x  8  matrix,  X  and  R  are  column  vectors  with  8  entries.  These  are  defined 
such  that  for  j  =  1,2, 


Xj  =  rHCi  - 


(173) 

(174) 


\z=:zii  for  k  =  1,2  , 


dz(k-i)  i 


(175) 


ana 


For  j  =  3, ...8, 


Gk^  =  e-'*''*"')  ^^1.=.,,  f/  for  k  =  3,4,  ..8. 

Xi  =  (Ci  - 

8- 

n.  .  =  _ I  f?  for  k  =  1  2 

'G’k,]  7  ®  lor  K  , 


(176) 

(177) 

(178) 

(179) 
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F;  for  k  =  3,4,..8. 


(180) 


Note  that  in  arriving  at  (172)-(180),  we  had  to  scale  the  original  continuity  equations  with 
appropriate  powers  of  7.  Further,  we  scaled  the  unknown  X  with  These  steps  are 

necessary  to  ensure  that  the  limiting  G  as  Rt  00  (i.e.  7  — +  00)  is  non-singular  and 
free  of  transcendental  terms  in  Rt. 

Denoting 

H  =  G-\  (181) 

it  is  clear  that  we  can  then  write  for  j  =  1,2, 


(182) 


while  for  j  =  3,. .8,  we  have 


C,  =  <7,  +  ^ 


(183) 


We  satisfy  (125),  (126),  (130)  and  (131)  at  z  —  2:0  and  (114)-(117)  at  z  =  Z2,  which  on  using 
(182)  and  (183),  is  equivalent  to  the  matrix  relation 


MZ^S, 


(184) 


where  A/  is  a  8  x  8  matrix,  while  Z  and  5  are  each  column  vectors  with  eight  entries.  These 
are  defined  such  that 

^  =  {Ci,C2A,C4A,Ce,C7,Csf,  (185) 


5=  (0,0,A„V3e  *'^0,0,-7a„^(z2),-7*c„p(2^2),-7'‘4^(^2),-7ftnp(2^2))^ 

+  {Pi,  P2,  P3,  Pa,  Psi  Pjy  Pa)^  , 

where  for  1  <  /  <  4, 

8  8  8 

P,  =  Y,  H,jRi  -  E  JWi.ie-"'*'"'’  E  HkjRi 

j=l  fe=6  j=l 

and  for  5  <  /  <  8, 

8  5  8 

j=l  k=3  >=1 

The  elements  of  the  matrix  M  can  be  written  ais 


(186) 


(187) 


(188) 


Mi,j  =  7  ^F-{zo),  for  j  =  1,2  and  Mij  =  Fj(zo)  for  j  =  3, 4..,  8,  (189) 
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^2,J  =7  forj  =  l,2  and  M2j  =  7  'Fj'(zo)  forj  =  3,4..,8 

=  +  +  forj  =  1,2, 

A^3.,  =  [7~‘A;‘F;“'  +  7’ViF;  +  y92A^fl,7-«FH  for  j  =  3, 4.., 8, 
=  [A;'F/  +  m,F/  +  m2F,']^  for  j  =  1,2, 


(190) 

(191) 

(192) 

(193) 


(195) 

(196) 

(197) 

(198) 


A/4,,  =  ^[Fj'  +  miF/l  +m2F/(2o)  for  j  =  3,4..,8  ,  (194) 

Ms^  =  F^(z2)  for  j  =  1,2,. .8,  (195) 

Mej  =  7~^F/(22),  for  j  =  1,2  and  Mej  =  Fj(z2)  for  j  =  3, 4.., 8,  (196) 

Mr,j  =  l-^Ff{z2),-  forj  =  l,2  and  M7j  =  7~’Ff(z2)  for  j  =  3,4.., 8,  (197) 

Msj  =  F,^(z2),  forj  =  l,2  and  Mgj  =  A^R,7’®FJ‘(z2)  for  j  =  3, 4..,  8.  (198) 

Note  that  in  arriving  at  the  specific  form  of  matrix  M,  which  is  nonsingular  in  the  asymptotic 
limit  of  Rt  —*  oo,  we  needed  to  multiply  each  of  the  original  boundary  conditions  by 
appropriate  scale  factors  involving  7.  Once  solution  Z  =  M~^S  is  found  to  (184),  the 
quantities  a„(2o)  and  6„(2o)  can  be  written  in  terms  of  Z  as: 

<!»(-ro)  =  7-'  {zi  +  ^  I 

+7-'  E  + 7-'  E  f;(^o)  (z,  +  e  .  (‘^s) 

i=3  i=6  \  k=i  / 

=  7''  |Pi‘(2o)Z,  +  Fi(zo)  {z2  +  7e-"'=l"''  E  «2.*ft)  } 

E  ^  E  ^(^0)  (^»  +  E  «,.»«»)  ■  (200) 

'7  2=3  7  j=6  \  fe=l  / 

In  terms  of  a„(2o)  and  6n(zo)j  as  shown  in  Appendix  A,  the  interfacial  deformation  coefficient 


(199) 


<i„  =  - 


ro'  +  mc^'Czo)  ’ 


and  the  radial  segregation  coefficient 


b„i{zo)  =  (^•C°'(2o)  -  cf  (Zo))  +  kbn{zo). 


(201) 


(202) 


Recall  that  the  function  reflects  the  common  exponential  part  of  the  growth  for 
the  functions  F“,  Fj  and  Ff.  As  discussed  earlier,  and  shown  explicitly  in  appendix  B,  for 
Rt  »  |/2c|,  with  Rt  »  1,  for  j  =3,4,5  decreases  exponentially  at  a  large  rate  with 
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increasing  z  from  a  value  of  1  at  z  =  zq.  This  is  uniformly  true  for  all  A„.  The  decreaise  is  at 
least  like  Further,  for  j  =  6,7,8  decreases  exponentially  with  decreasing  z 

from  a  value  of  1  at  z  =  22  at  a  rate  that  is  at  least  Thus,  the  elements  Mkj 

for  1  <  k  <  4,  6  <  j  <8  and  for  5  <  A:  <  8,  3  <  j  <5  are  which  is 

exponentially  small  in  Rt-  Further,  if  A„  >>  1,  it  is  clear  from  (157)  that  pi  ~  — A„  and 

P2  ~  A„  and  so  Mka  for  1  <  fc  <  4  and  Mk,\  for  5  <  A:  <  8  are  each  which  is 

transcendentally  small  in  A„.  Thus,  if  A„  or  z^  —  zq  is  of  the  order  of  some  positive  power  of 

Rt,  these  terms  will  also  be  transcendentally  small  in  Rt-  Further,  using  the  growing  and 
decreasing  properties  of  it  is  easy  to  see  that  the  the  elements  Pi  defined  in  (187),  (188) 
satisfy 

Pi  =  0  (203) 

Further,  the  term  X~^ appearing  in  the  third  element  of  S  in  (186)  is  transcen¬ 
dentally  small  in  A„  for  large  A„.  Noticing  the  special  structure  of  the  reduced  matrix  M, 
obtained  by  ignoring  transcendentally  small  elements  in  the  limit  of  large  A„  and  Rt,  it  is 
clear  that  the  solution  Z  to  (184)  has  components 


Zj  =  O  for  j  =  1,3, 4, 5, 

Further,  Z2  =  0(1)  and  each  Zj  for  j  =  6,7,8  are  at  best  0(1)  and  so 

eW,(zo)  =  O  , 

^  fijjR^  =  O  (e-2  for  ;  =  6, 7, 8. 


(204) 


(205) 


(206) 


Thus,  from  (199)-(206),  it  follows  that  the  effect  of  the  insulation  zone  is  to  exponentially 
quench  the  small  scale  components  of  the  raxlial  segregation  and  interfacial  deformation, 
which  correspond  to  b„i{zo)  and  d,,  for  large  An-  Most  of  the  contribution  in  the  summations 
(96)-(101)  therefore  comes  from  terms  where  A„  =  0(1). 

In  that  case,  the  matrices  G  and  M  simplify  further  in  the  eisymptotic  limit  Rt  >>  1. 
As  shown  in  appendix  C,  (184)  can  be  solved  in  closed  form  in  this  asymptotic  limit  resulting 
in  {Zq,  Zj,  Zs)  =  0(1), 


(Z3,  Z4,  Z5) 


1,  .V^  ,  I  .,/3' 
2'^’  6  ’  ’  2  ‘  6  , 


„-A„(zo-2/) 


(207) 


iZuZ-2) 


_ 2m2/73e-^"<^-^') _ 

(pi  -t-  mi  —  m2)  —  e(p>~P2h*2"*ol(p2  +  mi  —  m2)) 


(l,_ePi(-2-zo))  ,  (208) 
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K{zo)  ~ 


g-A„(20-2/) 


(r»'  +  mc"(zo)) 


(209) 


(210) 


(211) 


(212) 


where  the  effective  parameters  ^3,  mj  and  mj  are  as  defined  in  (134)-(136).  From  (199)  and 
(200)  and  the  above  asymptotic  results, 

an(-2o) - I  +  7 - ^ - 77-* - 7  (209) 

[pi  +  mi  —  m2  —  e(P‘“P2)l*2"*°l(p2  +  mi  —  m2)]J 

,  V  /  ^  mill  - 

K{zo)  ~  27-V3e-^'‘<"“-"'^7 - - 77 - 77— * - 7  (210) 

'  [pi  +  mi  —  m2  —  e(P»“P*)(**~*o)(p2  +  mi  —  m2)] 

From  (201),  we  find 

e-^n(2o-zi)  m2(l -m|){l 

dn  2/337  ^jiQ/  _|_  j^gO'^^io))  ^  ***  (pi  +  mi  —  m2  —  (p2  +  mi  —  m2)elP‘~P2h^2~^))  ^  ^ 

From  (202),  (210)  and  (211),  we  can  obtain  asymptotic  expression  for  6^(20)  as  well.  From 
(45),  (101)  and  recalling  that  Ji'(A„)  =  Jo(An),  we  get  interfaeial  slope 

/W  ~ 

7,  Ai'Vo{X„)  lt^‘{T«'  +  mc»'(jo)) 
m2(l -mf)(l 

X  1  +  7 - f-i - - ^77 - 7^  .  (212) 

(pi  +  mi  —  m2  —  (p2  +  mi  —  m2)e(P‘“P2)(**“*“l) 

Notice  that  in  putting  an  upper  limit  00  in  the  summation  in  (212),  we  have  to  go  through 
an  intermediate  analysis  where  we  replace  the  upper  limit  by  Ni,  where  Ni  »  1  but 
smaller  than  any  power  of  Rt  since  such  a  limitation  is  needed  for  the  validity  of  (159) 
used  in  deriving  (207)  and  (208).  The  contribution  to  the  summation  for  A„  even  larger 
is  transcendentally  small,  as  argued  earlier.  Now  since  Ni  is  far  larger  than  unity  and  the 
series  in  (212)  is  convergent,  the  leading  order  asymptotics  is  indeed  the  same  as  with  A^i 
replaced  by  00.  From  (44),  (100)  and  (202),  radial  segregation  in  the  solid  at  the  interface  is 

^(r  r  1  V'  4e^3(fec”'(^o)  -  C°'(^^o)) 

~  ;SAyVo(A„)  I^\T^'  +  m<P'{zo)) 

m2(l  -  m|)  {1  - 

^  (pi  +  mi  —  m2  —  (p2  +  mi  —  m2)e(p>~P2K**“^)) 

°°  JtiX  r)  4c/33m2fcae~'^"(*®“^^)  x  |l  — 

_j_  \  ^  / _ K _ j _  (213) 

n=i  Ay^Jo(A„)  (pi  +  mi  -  m2  -  (P2  +  mi  -  m2)e(Pi-P2)(^2-^)) 

Now,  consider  fluid  velocities  in  different  regions  of  the  melt.  By  using  (40),  (41),  (71), 
(72),  (98),  (108),  (150),  (154)-(156),  (182),  (183),  (185)  and  the  asymptotic  solution  (207)- 
(208)  and  simple  identities  of  Bessel  functions,  we  can  write  the  asymptotic  expressions  for 
radial  and  vertical  fluid  velocities  for  \z  —  zo]  =  0{Rj  '  )  as 


-{r,  zo) 


X  1-1- 


+  E 


(213) 


hi  xhMKir''  h 


(214) 
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(215) 


where 


ti  A^/Vo(A„)To'  ,t-3  ' 

^3  =  -l+i^  ,  N,  =  2,  Ns  =  -l-i^ 


3  '  ’  ■’  “  3 

This  gives  a  boundary  layer  a.t  z  =  zq  where  the  fluid  velocity  changes  rapidly.  Of  course,  the 
numbers  Nj  are  such  that  at  the  interface  z  =  zq,  each  of  these  expressions  reduce  to  «  =  0, 
iw  =  — 1.  Note  that  the  expressions  (214)  and  (215)  imply  that  the  velocities  change  rapidly 
within  a  boundary  layer  that  scales  as  Aside  from  the  factor  which  generally 


depends  on  n  and  the  heat  transfer  parameters,  it  is  easily  seen  that  within  the  boundary 


e  •^*1^  while 
g-Ai(*-2o)  The 


layer,  the  dimensionless  horizontal  velocity  u  scales  approximately  as 

the  vertical  nondimensional  fluid  velocity  lu  +  1  scales  approximately  as 
form  (214)  and  (215)  suggests  that  there  will  be  skinny  convection  cells  within  this  boundary 
layer.  More  detailed  discussions  on  the  parameter  dependence  of  velocities  will  be  taken  up 
in  section  7. 

Now,  let’s  consider  the  neighborhood  of  z  =  zjj,  the  end  of  the  insulation  zone  in  the 
melt.  We  first  consider  z  <  zjj.  Since  prior  analysis  shows  that  for  j  =  6,7,8,  each  of  Zj 
and  therefore  Cj  is  at  best  0(1),  it  follows  that  is  transcendentally  small.  Thus, 

from  (153),  (174),  (178),  (183)  and  (C15),  it  follows  that  for  j  =  6,7,8, 


k=l 


A„|-?,r°(z„)  +  ,,] 

KRt  +  A5 


(216) 


where  to  the  leading  order,  it  is  known  that  =  0(1).  Further,  since  Zj  for  j  =  3,4,5 
turned  out  0(1)  at  best,  it  follows  that  for  j  =  3,4,5,  will  be  transcendentally 

small.  Thus,  in  the  expression  (150),  for  the  stream  function  coeflBcient,  the  leading  order 
contribution  comes  only  from  j  =  6,7,8  when  z  <  zjj  with  \z  —  z//|  small  enough  to  be 
inside  the  boundary  layer  (actually  an  internal  layer).  It  is  clear  from  (150)  that 


j=6 


j=6 


Therefore,  from  (98)  and  (108), 
eif) 


7  /  \  \2A„«/?Te[  qiT^{zjj)  'f'  9^]  ij  'idz-zn) 


(217) 


n=l 


Here,  unlike  what  happens  near  z  =  zq,  most  of  the  contribution  to  the  summation  occurs 
when  A„  =  0{R^*)  »  1.  In  that  case  7,  determined  from  (158)  scales  as  and  hence 
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from  (71)  and  (72),  we  conclude  that  there  is  a  layer  that  extends  downwards  from  z  —  zu 
of  thickness  where  the  radial  velocity  u  scales  as 


eK/ii/*l-gtT°(z/j)  +  qi] 


JH)' 


(218) 


and  the  vertical  fluid  velocity,  u>  +  1  away  from  r  =  1  scales  as 

eK[-qiT°{zii)  +  92] 

J’O' 


(219) 


Near  the  corner,  where  each  of  1  —  r  and  z  —  2//  is  the  vertical  velocity  also  scales 

as  ‘n  (218),  i.e.  it  is  large  as  gravity  effects  become  large.  Similar  arguments  for  2  >  2// 
can  be  advanced  to  show  that  there  is  also  a  boundary  layer  that  extends  upwards 

of  2  =  2//  where  horizontal  velocities  scale  as  in  (218)  and  vertical  fluid  velocities  scale  as 
(219).  Within  this  layer  around  2  =  2//,  there  are  convection  rolls  that  can  be  deduced 
from  the  r  and  2  dependence  in  (217).  We  can  put  forward  arguments  near  2  =  22  as  well 
to  show  that  there  is  a  similar  boundary  layer  where  convection  occurs. 

Away  from  these  boundary  layers  at  2  =  20,  2  =  2//  and  2  =  22,  the  horizontal  and 
vertical  fluid  velocity  components  for  2  >  2//  scale  as  as  can  be  deduced  from  (98), 
(108),  (146)  and  (153).  This  stays  0(1)  as  with  increasing  gravity  effects.  However,  there  is 
in  this  case  also  a  boundary  layer  near  r  =  1  of  thickness  R^^^*  where  the  vertical  velocity 
scales  as  •  and  therefore  intensifies  with  increasing  gravity. 

For  2  <  2//,  but  away  from  boundary  layers  at  2  =  zjj  and  2  =  20,  it  follows  from  (150) 
and  the  behavior  (155),  (156)  of  FJ  for  j  =  3,  ...8  that  these  only  gives  transcendentally  small 
contribution  to  Cn.  The  only  sizable  contribution  comes  from  possibly  (7i  and  C2,  which  are 
determined  in  (208).  Thus,  it  follows  that  in  this  region,  the  fluid  velocity  scales  at  best  as 
and  is  reduced  by  increasing  gravity.  Unlike  the  core  region  for  2  >  zjj  there  is  no 
boundary  layer  in  this  case  at  the  side  walls  near  r  =  1 . 

In  order  for  our  leading  order  analysis  to  be  self  consistent,  it  is  not  necessary  that  scaled 
fluid  velocities  u  and  u;  +  1  be  much  smaller  than  unity.  Indeed,  we  obtain  the  same  leading 
order  result  if  the  advection  term  v  •  Vu  were  totally  dropped  in  (3).  In  order  for  us  to  be 
able  to  linearize  everywhere  in  the  fluid  field,  it  is  necessary  that  the  Reynolds  number  based 
on  the  largest  fluid  velocity  Rej  «  1.  Based  on  the  estimates  of  the  velocities  in  the 
boundary  layers  at  2  =  20  and  2  =  2//,  this  would  require  that  each  of 


l/JO' 

KR^*e[-qiT°{zii)  +  ^2] 
i/T®' 


(220) 

(221) 
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For  the  complete  validity  of  our  analysis,  e  needs  to  be  small  enough  so  that  both  conditions 
(220)  and  (221)  are  met.  In  addition,  we  need  that  the  interfacial  slope,  as  given  in  (212), 
is  all.  The  latter  does  not  limit  e  particularly  since  for  larger  gravity  levels,  as  on  earth, 
this  quantity  is  small  for  most  experimental  conditions  without  e  being  small. 

Equations  (212)  and  (213)  are  general  results  on  interfax:ial  deformation  and  radial  seg¬ 
regation  that  are  valid  whenever  |/?c|  <<  provided  the  constraints  listed  in  table 

(2)  axe  satisfied.  However,  this  result  is  complicated  by  the  dependence  through  effective 
parameters  03,  mi,  m2,  which  in  turn  are  complicated  functions  of  other  parameters.  It  is 
therefore  instructive  to  look  at  a  special  limit  that  is  applicable  to  many  experimental  set 
up. 

6.1  Special  case  of  Ds/zq  «  1  and  X\{z2  —  zq),  \\Zi,  Aizo  >>  1 

In  this  case,  as  mentioned  in  section  2,  c,-  ~  jC2,  c®’  and  its  derivatives  negligible  at 
z  =  Zo,c°'(zo)  ~  -(^  -  1)^,  c®"(-2o)  ~  (j  —  1)^.  Therefore,  from  (134)-(136), 


(1  —  fc)  m(l  —  k)c2 

^(1  -  k)c2 

aZ)2  (r®'-  i^(i-  1)) 

Thus,  in  this  case  the  interfacial  slope  expression  (212)  reduces  to 


f.  4Ji(A„r)^i 

[‘-rl 

f  217  +  \/4D»  +  V  ] 

• ,  (222) 

dr 

h  AyVo(A„) 

where 

X  _  _ 

«»e[?i.r®(2r/)  -  92,] 

(223) 

/4/®«[7’®'-(i-l)i^]’ 

^2  ■ 

lulf 

~  12.]' 

(224) 

Jf-  _ _ 

I}(1  -  *:)C2(1  - 

(225) 

"  aD2(r®'-(i-l)mc2/D)- 

The  expression  (213)  for  radial  segregation  in  the  crystal  at  the  interface  becomes 


dcs 
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(l-fe)a/ 

D  dr 
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4Ji(A„r)  ['  -  fel 

A,V^Jo(A„)  (^  -I-  +  X'l  -  ^  -  <^3) 

(226) 
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where 


S,  =  ^(1  -  fc)c2,  (227) 

Each  of  (222)  and  (226)  can  readily  be  integrated  with  respect  to  r  since  the  integral  of 
•^i(Anr)  is  —  Jo(A„r)/A„.  This  readily  gives  expression  for  the  interfacial  shape  /(r)  and 
concentration  in  the  solid  at  the  interface  c^^ZQ^r). 

We  now  turn  to  finding  simplified  expressions  for  the  velocities  near  the  boundary 

layer  at  ^  =  zq  for  the  simplified  case  of  this  subsection.  Using  (214),  (215)  and  simplified 
expression  for  0^,  we  obtain  that  each  mode,  i.e.  coefficient  of  Ji(A„r)  for  the  horizontal 
fluid  velocity  u  scales  as 


(228) 


The  coeflBcient  of  Jo(Anr)  in  such  a  series  representation  for  the  vertical  fluid  velocity  tn  +  1 
scales  as 


-  q,,\  [1  -  f]. 


(229) 


7  Discussion  of  Results  and  Conclusion 


This  is  a  problem  with  many  non-dimensional  parameters,  as  seen  in  Table  1.  It  is  interesting 
to  note  that  as  Rj  — »  oo,  the  leading  order  asymptotic  behavior  for  the  interfacial  slope 
and  radial  segregation  contain  far  fewer  parameters,  as  seen  in  (212)  and  (213).  However,  in 
applying  this  formulae,  it  should  be  pointed  out  that  the  next  order  correction  in  the  result 
is  0{RcIRt)  or  0(i2^^^^),  which  ever  is  larger. 

To  the  leading  order,  we  find  that  the  interfacial  slope  and  radial  segregation  scale  as 
Rj  '  ,  when  other  parameters  are  held  fixed.  Thus,  in  qualitative  agreement  with  numerical 
results  of  Chang  &  Brown^®,  we  find  that  moderately  large  gravity  is  worse  than  large  gravity 
when  it  comes  to  controlling  interfacial  slope  and  radial  segregation. 

To  avoid  a  very  detailed  discussion  of  all  the  physical  parameter  dependence  through  mi , 
m2  and  ^3,  we  now  restrict  ourselves  to  the  special  case  in  subsection  6.1.  It  appears  that 
most  experiments  satisfy  the  additional  restriction  placed  in  subsection  6.1  so  there  is  not 
much  loss  of  generality  in  doing  so. 

The  first  notable  observation  is  the  crucial  dependence  of  interfacial  properties  on  the 
quantity  Zq  —  zj.  This  is  so  in  the  general  case,  where  (212)  and  (213)  hold,  as  well  as  in 
the  special  case  where  (222),  (226),  (228)  and  (229)  hold.  Physically,  zq  —  z/  is  the  ratio  of 
the  distance  of  the  interface  from  the  end  of  the  insulation  zone  in  the  solid  to  the  radius 
of  the  cylinder.  Since  the  smallest  An  is  Aj  =  3.83,  it  is  clear  that  an  arrangement  that 
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would  make  Zo  —  zj  (note  zq  here  is  determined  from  (37)  when  St  a  moderately 

large  number,  say  2,  will  ensure  that  is  at  best  which  is  rather  small. 

Since  A2  =  7.016  and  other  A„  even  larger,  it  is  clear  that  under  most  such  experimental 
conditions,  the  dominant  contribution  to  the  series  (214),  (215),  (222)  and  (226)  will  come 
from  the  n  =  1  term.  The  profile  of  each  of  the  interface  /(r)  and  concentration  of  solute 
c,(r)  will  each  look  like  Jo(Air).  In  this  context,  it  is  interesting  to  note  that  we  found 
that  the  numerically  computed  Chang  h  Brown*®  interface  for  the  largest  thermal  Rayleigh 
number  quoted  in  their  paper  is  in  good  agreement  with  a  Jo(Air)  profile  even  when  their 
numerical  calculations  deal  with  the  fully  nonlinear  flow.  The  coeflBcients  however  could  not 
be  matched  as  their  boundary  conditions  at  r  =  1  is  diflFerent  from  ours.  Comparison  with 
the  Adornato-Brown*’’  profiles  is  not  as  favorable,  especially  near  r  =  1,  presumably  due  to 
the  absence  of  an  ampoule  in  our  analysis.  Returning  to  relations  (222),  (226),  it  is  clear 
that  if  most  of  the  contribution  in  the  series  comes  from  the  first  term,  there  will  be  an 
approximate  proportionality  between  the  radial  segration  ^  (r,  zq)  and  the  interfacial  slope 
/'(r)  with  the  coefficient  of  proportionality  equal  to 


D 


C-1 


^  (25  + 


(230) 


Earlier,  Coriell  h  Sekerka^^  hypothesized  such  a  relation  based  on  a  purely  diffusive  calcula¬ 
tion.  Brattkus  &  Davis^*,  without  any  insulation  zone,  find  no  necessary  relation.  We  find 
here  that  with  a  proper  insulation  zone  thickness,  there  is  an  approximate  proportionality 
between  the  two  though  the  constant  of  proportionality  in  (230)  is  different  from  Coriell  & 
Sererka^^  by  the  appearance  of  the  second  term  within  the  square  parentheses.  Since  Aj  is 
fairly  large,  if  the  Peclet  number  ^  is  not  large,  then  it  can  be  expected  that  in  an  approxi¬ 
mate  way,  the  second  term  within  the  square  parentheses  in  (230)  can  be  replaced  by  5^, 
which  can  again  be  small  for  many  materials.  This  may  explain  why  Coriell  &:  Sererka^^  got 
reasonable  agreement  with  some  experiment  data. 

Now,  let  us  discuss  the  physical  meaning  of  the  scale  parameter  that  appears  in  both 
(222)  and  (226).  First  note  from  (223)  that  becomes  large  when 


r»'  _  (1  _  l)mc,/C 


(231) 


approaches  zero.  However,  the  term  (231)  is  also  present  in  the  denominator  of  Sz  in  (225). 
Thus,  when  the  expression  (231)  approaches  zero,  both  (222)  and  (226)  approach  a  finite 
limit.  However,  the  expression  (231)  has  to  be  kept  positive  to  avoid  constitutional  super¬ 
cooling.  Indeed,  in  the  absence  of  capillarity,  the  condition  for  the  onset  of  Mullins-Sererka 
instability  of  the  basic  6  =  0  state  is  that  the  expression  in  (231)  is  zero.  From  the  boundary 
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condition  (13b),  the  quantity  (\q\,T^{zi)  —  ft,]  is  the  negative  of  the  raxlial  temperature 
gradient  in  the  solid  at  r  =  1,  2  =  z/.  Thus,  for  a  dilute  alloy,  the  ratio 

is  approximately  the  ratio  of  the  horizontal  temperature  gradients  in  the  solid  at  r  =  1, 
2  =  2/  to  the  vertical  temperature  grawlient  in  the  melt.  Thus,  =  rt{K,f  ,  where 

Ks/k  is  the  ratio  of  the  thermal  diffusivities  between  the  solid  and  the  melt. 

Now  consider  the  parameter  It  is  clearly  the  2  derivative  of  the  log  of  heat  loss  rate 
at  r  =  1,  2  =  2/.  If  dimensional  coordinate  2  were  used,  it  would  be  the  product  of  cylinder 
radius  a  and  the  logarithmic  derivative  with  respect  to  2  of  the  heat  loss  rate  evaluated 
at  the  edge  of  the  insulation  zone  in  the  solid  side.  This  is  clearly  a  property  that  can  be 
controlled  by  appropriate  design  of  the  heat  transfer  properties  of  the  Bridgman  apparatus. 
When  most  of  the  contribution  of  the  series  in  each  of  (222),  (226),  (228)  and  (229),  comes 
from  the  n  =  1  term,  as  is  the  case  with  an  appropriately  large  value  zq  —  2/,  then  the  special 
choice 

^2  =  Ai  =  3.83  (233) 

will  have  the  effect  of  minimizing  the  interfacial  slope,  radial  segregation  as  well  as  fluid 
velocities  in  the  boundary  layer  near  the  interface.  When  62  <  3.83,  the  interface  slope  in 
(222)  is  positive,  meaning  that  the  interface  will  then  be  bulged  towards  the  solid  as  reported 
in  previous  numerical  computations  (Chang  &  Brown'®  for  instance). 

The  role  of  the  parameter  S3  as  defined  in  (225),  is  more  complex.  For  sufficiently 
dilute  alloys,  clearly  since  C2  is  small,  S3  will  be  small.  In  that  case,  the  term  within  the 
curly  parentheses  in  (222)  reduces  to  unity  and  therefore  in  that  case  there  is  no  explicit 
dependence  of  the  interfacial  shape  on  the  segregation  coefficient  k  or  the  Peclet  number 
except  through  and  zq  —  zj. 

Even  in  a  general  case,  the  explicit  Peclet  number  l/D  dependence  in  (222)  and  (226)  is 
weak,  except  when  Peclet  number  is  comparable  or  larger  than  Aj  =  3.83 

We  now  discuss  fluid  motion.  As  discussed  earlier,  vigorous  fluid  motion  is  confined  to 
the  boundary  layers  near  2  =  20,  2:  =  2//  and  2  =  22  eis  well  as  near  the  side  walls  r  =  1  for 
2  >  2//.  Elsewhere,  the  motion  is  0(1)  in  the  bulk  for  2  >  zjj  and  o(l)  for  2  <  2/; 
as  gravity  is  increased.  This  qualitative  feature  appears  to  be  in  agreement  with  previous 
numerical  work. 

We  now  discuss  details  of  the  motion  within  the  '  boundary  layer  of  the  interface. 
Here  the  scale  of  the  horizontal  and  vertical  fluid  velocities  for  the  special  case  of  section 
6.1  are  given  by  (228)  and  (229).  Aside  from  the  factor  of  (1  —  S2l\n)i  whose  significance 
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has  been  already  discussed,  the  remaining  part  of  the  scale  factor,  when  written  in  more 
common  notation,  as  explained  in  Table  2  is  given  by 


—  Rt-^  Pr-'  r,  (234) 

/c 

for  the  horizontal  velocity  u.  For  the  vertical  fluid  velocity  u;  -|-  1 ,  the  scale  factor  is: 

Re-^  Pr-'  r,  (235) 

1C 

where  ^  is  the  ratio  of  solid  to  melt  diffusivities,  Re  is  the  Reynolds  number,  Pr  is  the 
Prandtl  number  and  the  parameter  r,  is  given  by 


-  92.] 


(236) 


Physically,  r,  is  the  ratio  of  the  horizontal  temperature  gradient  in  the  solid  at  r  =  1,  z  =  z/ 
to  the  vertical  temperature  gradient  in  the  melt  corresponding  to  the  Tiller  et  aP  solution. 
Note  r,  differs  from  rt  due  to  an  additional  term  in  the  denominator.  The  requirement  (220) 
for  dropping  nonlinear  terms  in  the  Navier-Stokes  equation  can  be  written  as 


—R^^  Pr-^  r,  «  1 


(237) 


However,  we  are  of  the  belief  that  (237)  is  too  stringent  since  within  the  boundary  layer, 
the  extra  z  derivative  in  the  Stokes  operator  part  of  the  Navier-Stokes  equation  introduces 
a  factor  that  scales  as  which  only  needs  to  be  larger  than  the  convective  terms,  whose 
relative  size  is  given  by  the  left  hand  side  of  (237).  This  will  have  to  be  confirmed  by  a 
nonlinear  analysis  in  the  future.  Nonetheless,  it  seems  that  through  a  safe  choice  of  (zq  —  z/), 
the  requirement  (236)  could  perhaps  be  met  in  experiment.  Of  course,  zq  —  z/  cannot  be 
chosen  arbitrarily  large  since  large  temperature  gradients  have  to  be  maintained  to  avoid 
supercooling. 

Now,  consider  fluid  motion  within  the  boundary  layer  around  z  =  z//.  The  scale  of  the 
horizontal  velocity  in  (218)  can  also  be  written  as 


where 


R^'*  Re-^  Pr- 


-  92] 


(238) 


(239) 


is  the  ratio  of  the  horizontal  temperature  gradient  at  r  =  1,  z  =  z//  to  the  Tiller  et  aP 
vertical  temperature  gradient  in  the  melt.  The  vertical  fluid  velocities  in  (219)  is  smaller 
than  the  horizontal  in  this  bounday  layer  by  a  factor  of  R^^,  except  near  r  =  1,  where  the 
fluid  velocity  components  are  of  comparable  magnitude.  The  requirement  (221)  for  dropping 


nonlinear  terms  may  be  harder  to  satisfy  in  an  experiment,  even  though  there  is  a  smaller 
power  of  Rj  in  the  scaling  compared  to  (237).  This  is  so  because  there  is  no  quenching  factor 
g-Ai(2o-*/)  However,  because  of  the  boundary  layer  structure  of  the  solution,  the  interaction 
of  the  interfacial  properties  at  z  =  Zq  with  fluid  motion  near  z  =  zn  or  z  =  Z2  is  weak.  This 
suggests  that  one  can  have  a  highly  nonlinear  flow  in  the  bulk  of  the  melt,  yet  if  the  fluid 
Reynold  number  near  z  =  zq  is  not  large,  i.e.  (237)  holds,  then  the  asymptotic  results  (222) 
and  (226)  will  also  hold.  In  that  case,  it  would  be  interesting  to  see  if  the  condition  (233) 
can  be  implemented  through  appropriate  design  of  the  Bridgman  apparatus. 

Again,  because  of  the  boundary  layer  structure  (i.e.  exponential  decay  of  the  modes 
away  from  the  boundary)  of  this  problem,  as  it  would  be  true  for  any  thermally  or  solu- 
tally  stabilized  configuration,  introduction  of  nonlinearity,  wheii  important,  is  expected  to 
be  mathematically  manageable  since  by  scaling  the  z  —  zohy  Rj  and  using  the  scale  infor¬ 
mation  on  velocities,  a  nonlinear  boundary  layer  equation  with  only  a  few  parameters  can  be 
formulated.  Through  appropriate  matching  with  the  linear  solution  in  the  nearly  stagnant 
core  of  the  insulated  zone,  a  global  solution  can  presumably  be  constructed.  This  will  be 
subject  of  further  investigation. 
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Appendix  A:  Derivation  of  effective  interfacial  boundary  conditions 

The  purpose  of  this  Appendix  is  to  to  derive  effective  boundary  conditions  on  the  melt 
variables  a.t  z  =  zq  by  using  (105),  (128),  (118)  and  (119)  together  with  the  interfacial 
conditions  (120)-(124). 

First,  recalling  that  7^  is  a  linear  function  as  in  (35),  it  is  clear  from  (128)  that  for 
0  <  z  <  Z], 

a^{z)  =  a„p(z)  +  (Al) 


where 


<*jnp(^)  —  .  (  "b  92.1- 


In  the  insulated  zone  zj  <  z  <  2o,  on  the  other  hand, 

asn{z)  =  (A3) 

The  continuity  of  a^n  and  its  first  derivative  ai  z  =  zj  implies  that  each  of  Bi  and  B2  can 
be  expressed  in  terms  of  B\  and  B2.  Further,  from  (118),  B2  can  be  solved  in  terms  of  B\. 
Going  through  the  algebra,  we  find  the  following  expression  for  a^n  in  the  insulated  zone 
zi  <  z  <  Zo". 

(^4) 


where 


—  2  0«np(^/)  +  ^ 


1  L  ^ 


Using  (A4)  for  a^(zo)  and  a',„(zo)  and  eliminating  Bi  between  them,  we  obtain 


Kni^o)  -  Xiasn{zo)  =  Fj 


where 


=A„ 


(1  + 


'•"(!  _e-2>‘nzo)  ^ 

,,  ,f-G'2  +  a,„p(0)e-''"*^ 

K.  =  - 1  _  m 

Similarly,  the  linear  homogeneous  equation  (105)  with  condition  (119)  can  be  solved  in  terms 
of  one  arbitrary  constant.  By  eliminating  this  constant  between  6',„(^o)  and  65,1(20),  we  find 


Kni^o)  -  X2bsn{zo)  =  0 


(AlO) 


where 

P3  — 

1  —  ’ 

(/111) 

where 

2D.  ‘ 

1/112) 

2D,  ^4D2'''^"‘ 

(■413) 

From  (120)-(122)  (with  do  =  0),  we  can  express 

asn(zo)  +  ^b^izo)  an(2o)  +  ^bnizo) 

T^'  +  m<^izo)  ~  ro'  +  mc®'(^o)  ‘ 

(AM) 

Equation  (120)  then  becomes 

anizo)  =  niia*n(^o)  +  nub^nizo), 

(■415) 

where 

Jf -T^ 

T^'+f<^'izoy 

(416) 

ma  [if  -  T®'] 

(417) 

From  (122)  and  (A14),  we  get 

bn{Zo)  —  Tl2ia»n(2^)  +  n22bsn{zo)y 

(418) 

where 

(419) 

1  m  [c"'(zo)  -  lcf(zo)] 

”22  =  I  +  r  1  • 

(A20) 

Thus, 

,  .  ”22an(^o)  -  nubn^Zo) 

^3n(Zo)  —  , 

”22”ll  ~  ”]2”21 

(421) 

J,  -”2ian(^o)  +  ”llfen(^o) 

”22”ll  ~  ”12”21 

(/122) 

Further,  from  the  simplification  of  (123)  into  a'^{zo)  —  ■^a'^{zo)  (due  to  linearity  of  T°  and 
T®)  and  (A7)  and  (A21),  it  follows  that 


-I-  02bn  =  (A23) 
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where 


^1  = 
^2  = 


XiK,n22 


A„«:(n22nii  —  ni2«2i) 
XiK,ni2 


A„K(n22>lll  —  «12>^2l) 

^  2A„k,  [-(?2  +  a^p(0)e“-'"*^  - 
^  K  [1  —  e~*^"*o] 

Using  expressions  for  wn,  ni2,  n2ij  «22i  Xi,  Gi  and  G2  in  the  above  and  carrying  out  the 
algebra,  we  get 


(>124) 

{A25) 

(A26) 


/9i  = 
02  = 


K,{Tf' +  mc<>'(zo)){\+e-^^’'^) 
K{mcP‘{zo)  +  T^')(l  —  ’ 

K,ma(T^'  -T^){\  + 


iA27) 

{A2S) 


03  = 


K0{mcf^(zQ)  +  r®')(l  — 

-  nqu-nM  -  -  ^7y(i + 

/cfi  —  ’ 

r429) 

Using  (124),  (A14),  (A22),  we  obtain  another  effective  boundary  condition  for  the  melt 
variables  of  the  form 

6'„(«o)  +  rnibn{zo)  +  m2an(zo)  =  0, 


where 


In  the  above. 


”^1  —  ^4  +  X20bi 
Tn2  =  0s  A  X207 


k-\  m - c»"(i«)) 
D  ^  (TO' +  mco'(zo)) 


0h  = 


D,n\i 


D(niin22  —  ni2n2i) 
a{T'-Vm,f{zo)) 

a 

07  — 


(A30) 

(A31) 

(A32) 

(A33) 

(A34) 

(A35) 

(A36) 


D{niin22  —  ni2W2i) 

Using  expressions  for  nn,  nu,  ”21  and  7122  in  the  above,  the  expressions  for  0s  and  07  simplify 
to: 

D,  (mc^ (zo)  + 

"  D(m<^(zo)  +  ro')  ’ 


(A37) 
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Appendix  B:  Derivation  of  the  asymptotic  form  of  approximate  general  solution 

The  purpose  of  this  Appendix  is  to  find  approximate  expression  for  c„p,  a„p,  6„p,  F',  Fj 
and  Fj  appearing  in  (146)-(152)  Rj  »  1  with  |/2c|  <<  Rt- 

It  is  possible  to  find  uniformly  valid  asymptotic  representation  of  the  solution  for  all  A„  in 
this  limit.  Nonetheless,  this  expression  is  too  complicated  and  not  very  suitable  for  algebraic 
manipulation  when  Rc  «  So,  we  first  derive  easier  expressions  for  the  asymptotic 

behavior  of  the  solution  for  Rt  »  1  and  A^|/lc|  <<  Clearly,  even  when 

l^cl  <<  /iy*,  this  condition  is  not  uniformly  valid  for  all  A„.  Thus,  for  large  enough  A„, 
we  have  to  use  the  more  general  expression;  nonetheless,  from  the  arguments  in  section  6, 
we  know  that  the  detailed  br  havior  of  the  solution  for  large  A„  is  not  necessary  to  conclude 
that  the  contribution  to  the  series  (96)-(101)  from  large  A„  is  indeed  negligible  and  so  the 
behavior  of  the  solution  for  (/2c |  <  /2r^*  can  be  deduced  from  the  results  in  the  following 
subsection 

Bl.  Subcase  Aj[|fic|  <<  7®»  Rt  »  1 

From  (143),  in  the  heat  zone  zn  <  z  <  zq,  we  have 

+  \lRcCxh{z)  -  \IRtC2  ]  c„  =  A  [-A„(-9,r®  +  92)] .  (Bl) 

To  the  leading  order,  we  expect 

Cn  -  cl  (B2) 

where  c°  satisfies 

[£,£2£3  -  >IRtC2  ]  c„  =  £2  [-An(-9ir®  +  92)]  .  (B3) 

Thus,  we  can  write 

[£,£3  -  AJRt  1  C„  =  +  ?j1  +  E  7-"  {£1  Nisb/l  -  AJ Hr}  >  (B4) 


i=i 

where  A\,  A2  are  arbitrary  constants  and 

7  =  A^H;/'  (B5) 

=  (B6) 

Hy]  =  D[y^  -  Ai:i  +  s  (67) 

£3(»1  =  (»"  -  Aj)"  +  liy^  -  \l)  (68) 
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I 


^  2D  ^4D2 


On  solving  {B4),  we  obtain, 


(fllO) 


4  =  4,  + 

j=l  j=3 

(Dll) 

WnA'TA 

.0  A„[-9iT"  +  d 

WAt  +  a;)  Wajat  +  asp’ 

and  7j  are  the  ix  independent  roots  of 

(D12) 

^3[7j]^i[7i]  ->IRt  =  0 

(D13) 

that  are  labelled  such  that  for  Rt  »  AJ[, 

7j  ~  It^i 

(D14) 

where  Uj  are  the  six  roots  of  unity  defined  by 

W3  =  — =  — wg 

(D15) 

W4  =  — 1  =  —Ur 

(D16) 

Wg  =  — =  —Us 

(D17) 

In  order  to  find  the  next  order  correction  in  the  asymptotic  expansion 

c»  c»  +  ci  +  - 

(D18) 

we  notice  that  satisfies 

A  [a  A  -  Ajflr]  4  =  -4«.A[*4,)-A;fl,7-’  E 

where 


j=i 


1 


i=3 


Pi -Pi 

1 

7i  =  7i-;p. 

Then  one  can  solve  (B 19)  in  a  standard  manner  to  find  that 


(D19) 

(D20) 

(B21) 


4  =  -  7  £ 


,ti^2[Pi]{Di[Pi]i3[Pi]-A^72r} 
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{B22) 


where  is  a  particular  solution  of  the  form 


Cnp  =  iM{z  -  zo)  +  Ai)h{z) 


that  satisfies 


c,  [£,£3  -  xiftr]  4,  =  -a;«.£,[a4,i. 

When  Rt  »  AJ[,  the  above  asymptotes  to 


1  ,  X  fi.  I,  ,  r>-„  0  Sc,  ill-f'lCl 

~  f  Kf  D<^n,+  I 

Thus,  if  we  write 

<=.  =  '^  +  'tl-'CiFf  +  j2CiFf 

j=i  i=3 


(B23) 


(B24) 


as  in  (146),  then 

~  c°p  +  c*p  +  ..  . 

Note  that  in  the  special  case  jRt  <<  A)[,  c,»p  and  its  derivatives  are  0{R^^),  which  follows 
from  (Bl2)  and  (B23).  This  scaling  property  is  used  in  section  6.  Further,  in  (B24), 

0_ 

KRt] 


~  I  {,  -  +  ..  ^  tor  i  =  1,2  ,  (B25) 


^  ~  { 1  ASflrl  +  ■■}  tor  J  -  3, 4, 5  ,  (B26) 

where  we  choose  new  arbitrary  constants 

Cl  =  7^1,  Ci  =  Cj  =  7®44y  for  j  =  3,4,5  Cj  =  for  j  =  6,7,8 

(B28) 

Using  (145a),  we  can  solve  for  a„  once  c„  is  known.  Such  a  solution  also  needs  to  be  consistent 
with  (144a).  We  find  that 

8 

(B29) 


where 


>=i 


®np  ~  ®np  A’  0„p  "I"  ..•  , 


(B30) 
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where 


0  A;|-<ir  +  ft]  Rt 

A;fiT  +  A$  .-lAJiiT  +  AJP  ’ 


and  a^p  is  a  particular  solution  of  the  form 


^hp  =  -  ^q)  +  Bi]  h{z) 


that  satisfies 


^i«np  =  . 

When  Rt  »  the  expression  simplifies  so  that 


,1 


KRMz) 


■^l[  B  Jc„  ,  r  »  .iO„ 


L2[-D-^] 


-11  ”p 


(B31) 

(532) 

(533) 

(534) 


Further, 

F?  ~  ePj(*-*(ai-a)) 
1 


,  _  fl.PiA(z)i„fel  (-i.,felz,3fel  +  DXIRt)  .  \ 

c/eTfti,[py)(i..[pili3N-A;«T)  7  ’ 


F«  ~  — -i - 

J  -Y?  — 

7; 


(535) 


5,p,h(;.)L,[7i] (--^.[7i)^3[7i]  +  DKRt)  ■  \  o.. 

55r7ii-2[7il  (ii(7i)i^3l7,l  -  KRr)  ‘7  ^  ’ 


(536) 


f;  ~ 


7/-A^ 


5cP,M^)^i[7il  (--i>2[7i]^>3[7j]  +  DKRt)  ,  \  678 

55T7ii^2(7i](/^i(7i]i3[7il-AW 

(537) 


Using  bn  =  —CzCn  —  a„,  we  can  write 


A.  =  ‘^  +  E7‘Q1^  +  A;ft7-'EC>A?  , 

j=l  3=3 


(538) 


where 


6np  =  -AfCnp]  -  Onp  •  (^39) 

Simplified  expressions  are  possible  for  Rt  »  AJ[,  which  we  do  not  care  to  write.  The  only 
property  that  will  be  of  importance  in  the  analysis  in  section  6  is  that  bnp  is  tnis  range  of  A,, 
is  0{Rcl RtiR^^)’,  as  is  a„p,  which  follows  from  (B30),  (B31)  and  (B34). 

Further,  in  the  general  case  in  (B38), 


>  D7»  7 

F*  +  „  for;  =  3, 4, 5  , 

^  7i^2[7i] 


for  j  =  1,2  ,  (540) 

(541) 


~  +  ..  for;  =6, 7, 8  (542) 

7;^2[7j] 

Now,  we  comment  on  the  validity  of  the  asymptotic  behavior  given  so  far.  The  behavior 
of  the  particular  solutions  a„p,  6„p  and  c„p  is  clearly  valid  anytime  \Rc\  «  Rj-  This  is 
uniformly  true  for  all  A„.  Further,  expressions  (B25),  (B35)  and  (B40)  is  also  consistent  in 
this  regime  for  all  A„.  However,  the  expressions  (B26)  and  (B27)  do  not  remain  valid  when 
is  the  same  order  or  larger  than  7®,  i.e.  This  is  because  from  definition, 

^il7.]i-3(7i]  =  =  7^ 

Thus,  for  large  Rj, 

^i[7jR3[7il  =  0(7®). 

In  order  that  the  second  term  in  (B26)  and  (B27)  be  smaller  than  the  first  term,  it  is  necessary 
that  |A^5c[  <<  7®.  This  is  the  origin  of  the  restriction  for  the  asymptotic  behavior  shown 
thus  far.  Note  that  the  condition  A^|5cl  7*  cannot  be  uniformly  valid  for  all  A„  even 
when  I /2c  I  <<  /2y®. 

B2.  General  form  of  Solution  for  |/2c|  «  Rt 

In  order  to  find  uniformly  valid  expression  for  all  A„  and  at  the  same  time  find  expressions 
valid  for  I /2c I  <<  Rt,  we  find  six  independent  solutions  to  the  associated  homogeneous 
equation  in  (Bl)  in  the  WKB  form.  These  solutions  will  replace  the  expressions  (B26)  and 
(B27).  Once  this  is  found  appropriate  expressions  can  easily  be  found  to  replace  (B36),(B37), 
(B41)  and  (B42),  which  are  also  invalid  in  the  general  case. 

We  consider  WKB  solution  of  the  form: 

=  for  j=3,4..,8  (543) 

We  will  think  of  the  relation  (B43)  as  exact  and  relate  the  WKB  approximate  behavior 
through  the  relation 

Wj  ~  Wo,  +  Wi,  (544) 

We  find  that  for  5r  »  1,  with  |5t  —  5c|  >>  1 ,  a  uniformly  valid  expression  for  Wb,  is 
given  by  the  following  expressions 

Wo,  =  -  +  \l]idz  (B45) 

tVo.  =  -  £l(^n«T  -  +  Xl\Uz  (B46) 

(B47) 
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VVo.  =  ly^lUr  -  +  Xl]idz 

Wo,  =  J\xlRr  -  Xl^Hz))^  +  X^Jdz 
Wo.  =  /jxlRr  -  Xl^h(z))ie‘^f  +  X^jijz 


(B48) 


(B49) 


(B50) 


lu  the  above,  we  proceed  with  the  understanding  that  the  principal  argument  is  being  used 
in  taking  the  squareroots  and  cube  roots.  Wi^  for  3  <  j  <  5  is  determined  by 


''i,  =  rKiz)dz 

JZQ  * 


and  for  6  <  j  <  8  ,  is  determined  by 


^1.  =  rn(z)dz 

Jzn 


(B51) 


(B52) 


where  each  W^.  are  determined  in  terms  of  the  corresponding  through 

WS  (27(H^  -  Ajf  +  24Ai;(<  -  AJ)"}  +  (^  +  j)  Wi(WS  -  Ajf  - 

m(wS  -  Xlf 

(B53) 

Further,  note  that  when  Rt  »  |/2c|  with  <<  Rj,  these  expressions  reduce  to 


(556) 


in  agreement  with  (B14).  Thus,  the  expression  (B43)  is  a  generalization  of  the  earlier 
formulae  for  for  j  =  3,4,  ..8 

Now  consider,  finding  Fj*  for  j  =  3, 4..8.  It  is  determined  as  a  particular  solution  to 


£,f;  =  -XIRtF^ 

Through  a  standard  dominant  balance  procedure,  we  obtsun  for  j  =  3,  ..8, 


ir« - 


'  {wf-xir 

Since  F*  =  —Fj  —  CsFJ,  it  follows  that  for  j  =  3,4,  ..8, 


(557) 


(558) 


f;  - 


wFTJi  [((^i  -  -  7*  -  2«  -  A^)^W^;  +  4WfWj\Wf  -  Xl)  +  -^(Wf  -  Xl 

^  (559) 


Notice  that  it  is  necessary  to  retain  all  the  terms  in  (B59)  even  when  the  first  two  terms 
within  the  parentheses  are  clearly  0(7®)  while  the  remaining  are  0(7^)  because  there  is 
cancellation  between  there  is  cancellation  these  leading  order  terms  because 

and  this  need  not  be  larger  than  7®.  Relations  (B43)  with  asymptotic  behavior  (B44)  and 
relations  (B45)-(B53)  replace  the  more  restricted  expression  (B25)  and  (B26)  when  is 
0(7®)  or  larger,  provided  Rj  »  \Rc\-  The  remaining  expressions  (B25),  (Br2)  remain 
valid.  As  far  as  o„  and  6„,  expressions  (B57)  and  (B58)  replace  (B36)-(B37)  and  (B41)-(B42) 
respectively,  other  expressions  (B35),  (B30),  (B31),  (B32),  (B39)  and  (B40)  still  remain  valid 
as  long  as  Rt  »  j/lcl- 
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Appendix  C:  Asymptotic  evaluation  of  solution  to  MZ  =  S 

The  purpose  of  this  Appendix  is  to  carry  out  the  asymptotic  evalilation  of  solution  to 
(184)  when  Rt  >>  1  and  A„  =  0(1).  In  this  case,  7  defined  by  (147)  is  >>  1  and  the 
aisymptotic  relation  (B159)  for  the  roots  7j  hold. 

First,  we  consider  the  elements  Pi  through  Pg  that  appear  in  the  expression  for  S  in 
(186).  Since  for  1  <  /  <  4,  the  elements  A//, *6"*^*^*^^^  scale  as  for  j  —  6,7,8 

and  is  therefore  transcendentally  small  in  Rt.  We  ignore  such  terms  in  Pi.  Similarly, 
for  5  <  /  <  8,  the  elements  scale  as  for  j  =  3,4,5  and  is 

transcendentally  small  in  Rt.  This  leaves  us  with  the  contribution  from  terms  multiplying 
Hij  and  in  (187)  amd  (188).  Thus,  it  is  necessairy  to  calculate  the  first  two  rows  of  the 
matrix  //,  which  is  the  inverse  of  G. 

First,  consider  the  simplification  of  the  matrix  G,  whose  elements  are  defined  by  (175), 
(176),  (179)  and  (180).  On  examination  of  (153)-(170),  it  is  clear  that 

G  ~  G“  +  O  .  (<^0 


where  for  Jb  =  1,2, 


„k-l 


=  Pi  for  j  =  1»2,  Ggj  =  -V^_  3  for  j  =  3,  ..8  , 

7j 


(G2) 


while  for  A:  =  3,  ..8, 


o(*-2)  /  X  (k-3) 

=  =  G2j=(7)  forj  =  3,..8. 


(C3) 


Now  consider  the  problem  of  determining  the  the  first  two  rows  of  =  CP~' .  This  can 
be  conveniently  done  by  finding  the  first  two  components  X[  and  X2  of  the  vector  X‘  that 
satisfies 

G»X'  =  £;'  =  {0,0,0,,,1,..0)’’,  (C4) 

the  only  nonzero  element  of  is  a  one  at  the  /-th  entry.  It  is  convenient  to  define  symbols 
Tfc  so  that 


isb)  =  (j'  -  Xir  +  -Xl)  =  '£  r.s* 


k=0 


Thus, 


1 

ro  =  ,  ri  = - 2.  ^  =  -2X1  ,  rg  =  -  ,  r4  =  1. 

V  V 


((75) 


((76) 


Also,  we  define  r*  =  0  for  negative  k.  Notice  that  (158)  implies  that 

1 


'l]-K 
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So,  when  the  first  two  components  of  the  vector  identity  (C4)  is  written  long  hand,  we  get 


=  i:  =  7-’  E  ’•*  L 

j=3  k=0  j=3 


(C7) 


8  4  8 

p,x[ + p,xi  =  Y.i'‘iMii\x‘i + 6,3  =  r'  E  +  ')^i + .  (cs) 

i=3  fc=0  >=3 

where  is  the  usual  Kronecker  delta  symbol.  From  the  {k  +  3)-rd  element  of  the  vector 
identity,  with  k  ranging  from  0  to  5,  it  follows  that 


j=3  i=l 


(C9) 


Using  this  in  (C7)  and  (C9),  we  get 


X‘+Xi^  EpiiafclJf'  +  7-‘r,-,  + 

Pi  A'j  +  Pa  A",  =  E  +  6,3, 

j=l 


from  which  it  follows  that 


Htj  =  X‘,  = 

-1 

II 

II 

-1  ■ 

~  ^1,2  ^_5  (P3^/-3  ~  n-4) 


P2  -Pi 


(P2  -  Pi) 


Pi ^1,1  -  ^  -s(Pin-3  -  n-4) 

Pi  -  P2  (Pl  -  P2) 


(CIO) 

(Cll) 

(C12) 

{C13) 


Now,  let’s  consider  the  other  elements  /f°,  for  k  >  Z.  Clearly,  from  (C9),  it  follows  that 
that  for  k  ranging  from  0  to  5, 


E'r-SyWi!,  =  -7-‘-'£|-'  +  6kJ-3 

i=3  j=l 


(C14) 


Using  the  results  (C12)  and  (C13),  it  is  clear  that  the  right  hand  side  of  (C14)  is  0(7“*)  for 
/  =  1,2  and  0(1)  for  /  =  3,4,  ...8.  Since  the  coefficient  of  each  term  on  the  left  hand  side  of 
(C14)  is  0(1)  in  the  asymptotic  limit  7  — »  00,  it  follows  that 


=  0(7“*)  for  j  =  1,2  and  Hfj  =  0(1)  for  j  =  3,  ...8  (CIS) 

Now,  since  from  (174),  ^**12),  (B30),  (B31),  (B34),  a„p  and  its  first  derivatives  are 
0{RcIRt,  while  Cnp  and  its  derivatives  are  0(/?^*),  it  follows  from  (174)  and  (178)  that 
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each  of  Ri  and  are  O^Rd RtiI  ®),  while  R^  and  R4  are  0(7  *),  and  0(7  respectively. 
All  other  components  are  0(7“^  or  smaller.  This  means  that  for  j  =  1,2, 


8 


Y;  Hj,kRk  =  x 

((716) 

k=\ 

Pt  =  X  0  ('y-^XlRc'r'^] 

)  for  /  =  1, 2,3,4  . 

((717) 

Pi  =  X  0  {r^xlRc^r") 

for  /  =  5,6, 7, 8  . 

((718) 

Further,  we  notice  that  in  (186)  in  the  terms  beside  P/,  all  except  are  0(7"®Pc) 7~^)  = 
0(1),  where  as  S3  contains  a  term  that  scales  as  which  will  dominate  every  other 

term  for  An  =  0(1)  and  zq  —  zj  =  0(1).  Thus,  to  the  leading  order 

Here  we  comment  that  even  if  we  were  not  to  ignore  the  contribution  to  Si  for  /  >  5 
on  the  grounds  that  zq  —  zi  is  large  and  therefore  small,  it  would  not  affect  our 

leading  order  results  for  Zi,  Z3  through  Z5  because  of  the  the  special  structure  of  the  reduced 
matrix  M  in  this  limit.  Therefore,  the  results  quoted  in  (209)  and  (210)  would  equally  be 
valid  in  this  case  since  the  contribution  of  Zj,  Z^  through  Zg  for  interfacial  properties  an(^o) 
and  6n(-2o)  (given  by  (199)  and  (200))  are  exponentially  small. 

Now,  consider  simplification  of  the  matrix  M,  whose  elements  are  shown  in  (189)-(198). 
We  notice  that 

M  ~  (C20) 

where  the  only  nonzero  elements  of  A/°  are 


for  k  =  1,2  j  =  3,4,5,  (C21) 

K  =  .  (C22) 

^4,t  =  “Pi  -  ”>1  + 1"!  ,  AfJj  =  (-P,  -  m,  +  m,)  ,  for  j  =  3, 4, 5  , 

(C23) 

<1  =  e"'”'"’ ,  «?.,=  !  ,  M“^  =  -ui-''forj  =  6,7,8,  {C24) 

Mej  =  1  and  Mtj  =  Wj  for  j  =  6, 7, 8  ,  ((725) 

<1  = =  (C26) 

The  solution  to 

A/®Z°  =  5®  (C27) 
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can  be  found  in  closed  form.  We  find  (Zg,  Z®,  Z®)  =  (0,0,0), 


^0  7O  '7Q\  _  (  ^ 


4 -f) 


(028) 


Once  this  is  found,  it  is  clear  that 


{Z®,Z®)  =  - 


_ _ 

(pi  +  mi  —  m2)  —  e(p‘~P2)(**~^)(p2  +  mi  —  m2)) 


„Pl{z2-Zo)^ 


(C29) 
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Figure  and  Table  Captions 
Figure  1:  Schematic  diagram  of  the  Bridgman  apparatus. 

Table  1:  Nomenclature.  Note  all  listed  variables  are  non-dimensionalized. 

Table  2:  Relation  of  non-standard  symbols  with  more  common  notation  in  Literature. 
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Table  1:  Nomenclature 


T 

Nondimensional  temperature,  where  T:  temp.;  To:  planar  interfacial  temp  for  pure  material 

and  T}*.  top  temperature. 

c 

relative  concentration  of  one  binary  component 

V  =  (u,  w) 

Fluid  velocity  relative  to  velocity  of  container,  u  and  w:  radial  and  axial  components 

z 

,  where  z:  axial  coordinate  and  zi :  z  at  the  lower  cylinder  end. 

Xt 

Value  of  z  at  the  lower  end  of  the  insulation  zone  in  Fig.  1 

2II 

Value  of  z  at  the  upper  end  of  the  insulation  zone  in  Fig.  1 

20 

Value  of  z  at  the  unperturbed  interface 

Value  of  z  at  the  cylinder  top  in  Fig.  1 

r 

2  ,  where  f:  radial  coordinate. 

a 

Nondimensional  thermal  expansion:  a  =  d(7'3  —  To),  where  a:  volumetric  thermal  expansion  coefficient. 

Coefficient  of  volumetric  expansion  due  to  increase  in  c 

u 

Actual  Visco8ity/(at/a)  ,  where  a:  cylinder  radius,  Ua  container  velocity. 

K 

Thermal  diffusivity  in  the  melt  /{alia) 

D 

sclutal  diffusivity  in  the  melt  /(aUa) 

K, 

Thermal  diffusivity  in  the  solid  /{aUa) 

D, 

solutal  diffusivity  in  the  solid  /{aU^) 

Biot  number  related  to  horizontal  heat  transfer  on  the  melt  side 

,  where  (  «  1  is  chosen  so  that  91  =  0(1) 

Biot  number  related  to  horizontal  heat  transfer  on  the  solid  side 

,  where  (  «  1  is  chosen  so  that  91,  =  0(1) 

m 

Nusselt  number  related  to  horizontal  heat  transfer  on  the  melt  side 

f?2. 

Nusselt  number  related  to  horizontal  heat  transfer  on  the  solid  side 

C2 

Concentration  at  the  cylinder  top 

k 

Seggregation  coefficient 

p 

,  where  po:  density  at  temperature  To  and  concentration  C2 

9 

acceleration  due  to  gravity  x  a/Ug 

Si 

Stefan  number: 

[T2—To)C, 

do 

surface  tension/(Latent  heat  x  a) 

m 

Slope  of  liquidus  line  x  (T2  —  To) 

Rt 

Thermal  Rayleigh  Number 

where  T®  :  z  derivative  of  T°,  the  unperturbed  T  corresponding  to  <  =  0 

Rc 

Solutal  Rayleigh  Number 

where  c°:  interfacial  concentration  of  e  =  0  state 

Kz) 

Convenient  symbol  for  e“ 

A„ 

nth  positive  zero  of  Bessel  function  Ji 

7 

Convenient  symbol  for 

Table  1:  Nomenclature.  Note  all  listed  variables  are  non-dimensionalized. 
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Table  2:  Relation  of  non-standard  with  standard  parameters  in 

Literature 


D 

=  Pc“L  where  Pe\  Peclet  number 

D-'  «  . 

D. 

=  {D,/D)Pe~^ ,  where  D,/D:  Ratio  of  solid  to  melt  solute  diifusivity 

No  restriction  . 

1/ 

=  where  Re:  Reynolds  number  based  on  translation  rate  Ua  ■ 

Re  «  . 

Zj 

=  A~^,  where  A:  Ratio  between  cylinder  radius  and  length. 

No  restriction. 

1C 

=  Re~^ Pr~^ ,  where  Pr  =  Prandtl  number 

/c/22  >>  1- 

K, 

=  (/c,/K)Rc~*Pr“*,  where  k,/k:  Ratio  of  solid  and  melt  heat  diffusivity. 

'e./z2  »  1- 

Rt 

Thermal  Rayleigh  Number  ,  where'denotes  dimensional  quantities, 

T®:  Limiting  linear  Temp,  profile  of  Tiller  et  al  (1953)  for  /c/22  »  1- 

Rt  »  1. 

Rc 

Solutal  Rayleigh  number  ^bere  c^:  Tiller  et  al  (1953)  interface  concentration. 

|Rc|  « 

ga 

=  RTRe~^Pr~^A-\  where  A  = 

No  expl.  restr. 

gP 

=  RTRt-^Pe-\ci  -  C2)-‘A-‘,  a  =  1  - 

No  expl.  restr. 

Table  2:  Relation  of  non-standard  symbols  with  more  common  notation  in  Literature. 


53 


REPORT  DOCUMENTATION  PAGE 


Form  Approntod 
Om  No.  OAM-OfM 


rulilli  iwnninfl  tiiiriin  ^n^1^ll  nnttinttir  irf  '~*i — i' —  i  r*~ — T~‘  —  Mr*9»  i  Hour  ocf  r«won*.  mciudiis  Ihc  tunc  (or  mwouto  rnttruciiom.  scorrhuif  ciisiing  dou  soutrn. 
oMlwiiia  «Mf  fTumrimtitii  llic  dM  and  eowpittiiig  and  rawtawing  th«  <o<t«<ttoA  of  mfonnaooA  tend  coffMoanto  raaording  im  burden  estimate  or  any  other  aspect  of  ttus 

Sieepon  of  inform  mow.  Niduding  luggeifioni  for  rodurtwn  ctus  burden,  to  Washings  Headduarten  Services,  drectoraie  for  wtormation  Operations  and  Reports.  I2  >S  ieffenon 
Ooem  lOiAeiby.  Sswfd  130*.  Artiwgfon.  VA2a202’43M.  and  to  the  Office  of  Management  and  tudget.  eapemvorfc  deduction  Progea  <0704-01 W.  drashmgton.  DC  20S03 


1.  ACINCV  USf  ONLY  Oomit  bUtik)  1 2.  lUrOIIT  OATE  1 3.  AEKMT  TYPE  AM)  OATES  COVEMO 

I  September  1993  I  Contractor  Report 


a  TITU  AND  SIWTITU  S.  EUNOWG  NUMKKS 

CONVECTION  EFFECTS  (»l  RADIAL  SEGREGATION  AND  CRYSTAL  C  NAS1-I860S 

MELT  INTERFACE  IN  VERTICAL  BRIDOIAN  GROWTH  C  NASl-19480 


7.  raVOIIMMG  ONGAMZATION  NAME(S)  AND  AOOMSS(ES) 
Institute  for  Computer  Applications  In  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 
Hatton,  VA  23681-0001 


S.  SPONSOSSIG/MONnOlimG  A6ENCV  MAME(S)  AND  AOOAESS(ES) 
National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  VA  23681-0001 


WU  505-90-52-01 


S.  KRFOAMWG  ORGANUATION 
AEPOST  NUMSEK 


ICASE  Report  No.  93-72 


10.  spoNSOamc  /  monitoring 

AGENCY  REPORT  NUMSER 

NASA  CR-191541 
ICASE  Report  Mo.  93-72 


11.  SUPPLEMENTARY  NOTES 
Langley  Technical  Monitor: 
Final  Report 


Michael  F.  Card 


12a.  OBTRMUnON/AVAAASNJTY  STATEMENT 

Onclasslfied  -  Unlimited 


Submitted  ti  Physics 
of  Fluids  A 


12b.  DISTRIBUTION  CODE 


Subject  Category  34 


13.  abstract  CMaxini«Nn200«»ardU 

He  analytically  study  the  Influence  of  convection  caused  by  horizontal  heat  transfer 
through  the  sides  of  a  vertical  Bridgman  apparatus.  We  consider  the  case  when  the 
heat  transfer  across  the  side  walls  Is  snail  so  that  the  resulting  Interfaclal  de¬ 
formation  and  fluid  velocities  are  also  small.  This  allows  us  to  linearize  the 
Navler-Scokes  equations  and  express  the  interfaclal  conditions  about  a  planar  Inter¬ 
face  through  a  Taylor  expansion.  Using  a  no  tangential  stress  conditions  on  the 
side  walls,  asymptotic  expressions  for  both  the  interfaclal  slope  and  radial  segreg¬ 
ation  at  the  crystal-melt  Interface  are  obtained  In  closed  form  In  the  limit  of 
large  thermal  Rayleigh  number.  It  Is  suggested  that  these  can  be  reduced  by  appro¬ 
priately  controlling  a  specific  heat  transfer  property  at  the  edge  of  the  Insulation 
zone  in  the  solid  side. 


IS.  NUMBER  OF  PAGES 

58 


14.  SUBJECT  TERMS 

crystal  growth;  convection  effects;  Bridgman  apparatus 


17.  SECURITV  CLASSIFICATION  IS.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSffKATION  20.  LIMITATION  OF  ABSTRAa 
OP  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Thiclssslfled 


NSN  7S40-01-280-5500 


'4'IJS.  GOVERNMENT  PRINTING  OmCE:  l««3  ■  52S4«4/M07S 


Standard  Form  298  (Rev  2.89) 

Pnwnbcd  b,  ANSI  Std  Z»-<t 

n«-(02 


